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nroberties  oi  the  cloud  of  such  fiber  particles. 
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FOREWORD 


Th  e  formal  purpose  of  this  report  is  to  address  the 
Statement  of  Work  of  Contract  No.  F496 20-84-C-Q045  .  Considerable 
additional  effort  has  been  expended  in  Section  I  to  provide  a 
comprehensive  presentation  of  onr  most  recent  theoretical  results 
in  electromagnetic  scattering  and  absorption.  Since  this  work 
enoompasees  over  20  years  of  effort  (see  for  example,  N.E. 
Pedersen,  J.C.  Pedersen,  H.A.  Bathe,  "A  New  Method  of  Radar 
Target  Concealment,"  Proc.  Tri-Services  Radar  Symposium,  San 
Diego,  1969),  it  has  not  been  possible  to  include  many  topics, 
such  as  comparison  of  baokscatter  predictions  with  experimental 
re  s  ul  t  a . 

In  Section  II  we  address  the  question  of  eleotit? agneti c 
absorption  and  subsequent  thermal  radiation  by  very  small 
partiolea.  As  in  Section  I,  we  have  attempted  to  provide  a 
o on prehen ; ive  treatment  of  this  subject.  The  material  presented, 
together  with  that  in  the  cited  references,  permits  the  solution 
of  a  wide  range  of  related  problems. 

During  the  course  of  the  program,  a  Hewlett  Packard  Model 
9000/520  computer  was  provided  by  AFOSR.  Without  the  use  of  this 
fast  machine,  a  large  fraction  of  the  computations  presented 
herein  would  not  have  been  possible.  It  is  presently  being  used 
in  conjunction  with  other  DoD  programs,  including  AFOSR.  We  are 
very  grateful  for  this  contribution  by  AFOSR,  and  believe  that  it 
has  contributed  greatly  to  DoD  and  to  our  ultimate  understanding 
of  aany  of  the  phenomena  with  which  we  have  been  dealing. 
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SBCTION  I.  BLBCTB0HA8HBTIC  TBBOBT 


1.1  Ivtradaati 


Tki  groundwork  for  the  problea  of  ioattirii|  of  eleetro- 
aagaetie  vivm  by  thla  wirti  wit  laid  by  Poo  11  imgt  oa2 ,  aad  lotor 
Holloa2,  *ke  developed  tho  governing  latogral  oqaotioaa.  la 
1947,  tho  olaaaie  papo r  of  Vaa  Vlaek,  Block  aad  Haaeraesh 
proaoatod  extensive  roaalta  for  ba o k a e a t to r 1 ag  by  porfootly 
ooadaotiag  virea,  aaiag  Halloa's  equation.3  Shortly  aftor  that, 
a  variational  aathod  baaod  oa  Po o k 1 i a « t o a • a  equation  vaa 
doaoribod  by  Tai*,  yielding  roaalta  for  tko  aoat  part  la 
agrooaoat  with  tkoao  of  Vaa  Vloek  at  al.  Tko  wariatioaal 
technique  vaa  lator  oxtoadod  to  viroa  of  flaito  conductivity  by 
Caaaody  aad  Palaborg.  althoagk  tkoy  eoatiaaod  to  coaaidor  oaly 
baokaoa ttoriag  at  aorval  incidence.3 

la  rolatod  vork,  a  aaao r i c a  1 1 y- or 1 o at o d  aolutioa  of  tko 
Pookliagtoa  equation  vaa  oarriod  oat  by  Kiehaoad  for  both 
porfootly  conducting^  aad  f  iai ti-oonduct iy ity2  viroa.  A  aiailar 
approaoh  vaa  takoa  by  Medgy oai-Mi t sehaag  aad  Bftiaiu.*  Applic¬ 
ability  of  a  Vioaor-Hopf  toehaiqao  vaa  akota  by  Chen. 9  Aa 
iatoroatiag  approaoh  ia  toras  of  oatgoiag  aad  reflected  viroa  haa 
booa  givoa  by  Shea,  Va  aad  King,  althoagk  oaly  Halted  roaalta 
are  available.10  Curved  vires  have  also  booa  eoasidored,  tho 
ciroalar  loop  by  Ko ay o aa j i aa2 1  tad  tho  aoro  goaoral  naao  by 
Mei.22  A  good  reviev  of  tko  litorataro  ia  givoa  by  Einarssoa.23 

Tho  proaoat  authors  have  otadiod  this  probloa  ia  a  series  of 


>••  tki  vttiitioul  attkod,  to  that  the  work  eta  bt  regarded  tt 
am  eattatioa  of  tht  work  of  Tai*  tad  Caaaedy  tad  FaiabergS  to 
lieltda  botk  fialte  conductivity  tad  arbitrary  tagle  of 
iaeideaee.  fioiifir,  wt  f  iad  it  ton  convenient  to  derive  tkt 
bttie  equations  ia  tkt  eoataat  of  Otltrkia't  attkod,  wkiek  Joata 
kta  tkowa  to  bo  exactly  equivalent.** 

Diffortatial  aeattoriag  patteraa  art  tkoa  eoapated,  tloag 
witk  tkt  aeattoriag,  abaorptioa  tad  oxtiaotioa.  Theta  Croat 
aeetioaa  are  obttiaed  by  iategratiag  tkt  aoratl  eoapoaaat  of  tkt 
Poyating  veotor  over  tka  aarftca  of  tka  fibar,  eatbliag  at  to 
avoid  tka  iatagrttioa  ovat  tka  ftr-field  apkara  ttatlly  aaployad 
to  eoapate  aeattoriag.  Ia  additioa,  wa  fiad  that  eaargy 
ooa t i da r a t i oa t  art  axaetly  aatiafiad:  tka  oxtiaotioa  eroat 

aaotloa,  wkiek  by  the  optieal  tkeoraa  watt  equal  tke  laogiaary 
part  of  tka  forward  aaplitada,  it  idaatioally  equal  to  tka  tu  of 
tke  abaorptioa  tad  aeattoriag  erott  aeetioaa.  Tkla  raaalt  it 
partiealarly  rateable  for  thoaa  applicatioaa  ia  wkiek  wa  study 
tka  aeattoriag  and  abaorptioa  propartiaa  of  a  cloud  of  a  uek  fibar 
par  tielaa. 

Baoauaa  of  tka  approxiaata  nature  of  the  variatioaal 
approach,  it  ia  iaportant  to  baek  up  the  reunite  by  independent 
ooaputatioaa  whereror  poaaiblo,  in  addition  to  ooapariaon  with 
existing  raaulta.  In  tka  Rayleigh  region,  where  fibar  length  ia 
aaall  ooaparad  to  waralaagth,  el  a  e  t  r  o  a  t  a  ti  c  oo  aai  da  ra  ti  out  would 
be  expected  to  doainata.  A  quaaistatic  nodal  ia  dereloped  and 
ekaekad  againat  variational  reaulta.  At  high  frequenciea,  on  the 
other  hand,  tka  abaorptioa  croaa  aection  par  unit  length  auit 


Approach  that  cf  the  infinitely  long  fiber.  This  latter  quantity 
i»  ooapo';sd,  and  agreement  is  seen  to  be  good. 

This  work  was  originally  begun  for  applications  in  the 
microwave  region.  In  order  that  the  computations  nay  be  extended 
into  the  infrared  and  visible  regimes,  it  is  vital  to  incorporate 
the  optical  properties  of  the  fibers.  We  do  this  for  applicable 
materials  by  employing  the  Drude  model  for  conductivity  (or 
complex  dielectric  constant),  and  include  at  well  the  dependence 
of  conductivity  on  both  fiber  diameter  and  electron  mean  free 
pa  th. 


1.2  Theory 


Consider  the  thin  conductive  fiber  of  radius  a  and  half- 
length  h,  as  shown  in  Fig.  1-1,  hazing  length  to  diameter  ratio 
h/a  >>  1.  An  electromagnetic  wave  is  incident  along  a  direction 
making  an  angle  6^  with  the  positive  z-axis,  as  shown,  and  we 
suppose  that  the  fiber  radius  is  very  small  compared  to  incident 
wavelength,  i.e.  ka  <<  1,  so  that  the  resultant  scattered  wave 

has  rotational  symmetry. 


With  t  i  m  e- de  p  e  nd  e  n  c  e  factor  exp(-iut)  suppressed,  the 

12  4  5 

problem  is  described  by  the  integral  equation  '  *  * 


X+  h 

dx*  I(x')  K(z* ,x) .  (1-1) 

h 


The  left  hand  side  of  this  equation  expresses  the  axial  component 

of  the  E-vector  along  the  surface  of  the  fiber  in  terms  of  an 

induced  line  current  I  ( z)  (equal  to  2na  times  the  induced 

2  5  15 

surface  current)  and  a  surface  impedance  '  ' 


who  r  e 


+  io>n 

1,8  =  urn 

K  =  (  0)^  |i  8 


JQ  (ta) 
Jl  (xa) 

iupo)  1  ^ 


(1-2) 


(1-3) 


is  the  compiwi'  propagation  constant  within  the  fiber  (we  will 
assume  free-space  permeability,  however)  and  Jq,  J  ^  are  Bessel 
functions.  It  is  also  tacitly  assumed  that  lx/kl^  >>  1,  so  that 

4 


radial  variations  of  the  field*  within  the  fiber  are  much  aore 
rapid  than  axial  ones.  The  first  tern  on  the  right  hand  side  of 
Eq.  (1-1)  in  simply  the  axial  component  of  the  incident  E  field. 
Finally,  the  integral  tern  gives  the  E  field  contribution  due  to 
the  induced  current,  where  the  kernel  is  given  by 

I  (z'.x)  -  (1  +  d2/k2dz2)  ( 1/ 2n )  ]  dO'  (1/R)  eikR,  (l-4a) 

0 

r2  a  (  x—  z  *  )  2  +  4  a2  sin2  ( 0  *  / 2  >  .  (l-4b) 


The  ourrent  is  now  written  as  a  linear  coabination  of  the 
trial  functions 


e  ^ikz  .  e  ±-ikz  co*  ®i.  (1-5) 

Note  that  the  first  two  of  these  represent  resonant  currents  that 
can  exist  on  the  infinite,  perfectly  conducting  fiber  with  no 
external  field  present*,  and  should  continue  to  be  appropriate 
for  moderate  to  large  conductivity.  The  term  exp  (  ikz  cos  0  ^ ) 
gives  precisely  the  forced  response  of  the  infinite  fiber 
regardless  of  the  value  of  conductivity.  Finally,  the  term 
exp  (-ikz  cos  0*)  is  included  so  that  the  trial  functions  can 
satisfy  the  boundary  condition 

I  (+  h)  »  0  (1-6) 
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of  vanishing  current  a c  the  and*  of  the  fibor  tad  still  Maintain 
sons  flexibility. 

Following  this  prescription,  the  surface  current  takes  the 

fora 

IU)  »  ae  f  0  ( x )  +  a0  f0(x) 

“  a,  [ooc  kx  cos  (kh  cos  8^)  -  cos  kh  cos  (kx  oos  6j)] 

+  a0  Iain  kx  sin  (kh  cos  8^)  -  sir  kh  sin  (kx  cos  Ojil.  (1-7) 


In  order  to  deteraine  the  coefficients  by  the  Oalerkin 
aethod,  first  substitute  Bq.  (1-7)  into  Eq.  (1-1),  then  Multiply 
separately  by  fe(x)  and  f0(x)  and  integrate  over  fiber  length. 
Fron  the  two  resulting  equations  one  finds  that 


4niE 


o  “  e ,  o  (even,  odd)  (1-8) 


in  terms  of  the  ( dinen s io nl e s s )  quantities 


g0  -  k  sin  8 i  /  dx  fc 

-h 


r* f„<,>  "* 

•'-it 


+  h  » 

■  -4ni  (q»/q0)  f  dt  fc(x) 


r 

-  h 


t 


ds'fa(s')  K*'.*) 


(1-9) 


/•+k  /•►h 

Ta  “  k  I  d*  fo<  *>  I 

*'-h  ,'-h 


for  a  ■  •,  o. 


All  of  these  Integra’  i  can  b#  evaluated  analytically  under  the 
approxins  tion  ka  <<  1;  the  results  are  ratber  lengthy  and  *e  do 
not  give  then  here.  Note  that  sons  of  these  expressions  were 
evaluated  by  Tai,  for  a  soaevhat  simpler  kernel.*'  **  As  noted 
earlier,  tho  variational  procedure  used  by  Taid ,  as  well  as 
Cassedy  and  Faiabergd,  would  give  identical  results. 

The  scattered  field  can  now  be  expressed  ns  an  integral  of 
the  induced  current.  For  far-field  scattering  in  the  O-direction 
(see  Fig.  1-1)  one  ka*^ 


B,  -  0  Bi  (1/kr )  eikr  S(0i,0), 


(1-10) 


with  far-field  amplitude  given  by 


+  h 

S  (0i,  0)  -  (kn0/Ei)  sin  0  |  dz  I(z)  e~ ikz  co#®.  (1-11) 


/ 

-  h 


For  the  differential  (or  bistatic)  oross  section  one  has 

od  *  (4n/k2)  \  S  (Oi.0)  I2,  (1-12) 

and  similarly  the  bac  kc  ca  t  ter  ing  or  radbr  cross  section  is  given 
by  (again  see  Fig.  1-1) 
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“  (4n/k2)  I  S(8i,*-0£)  I2. 


(1-13) 


Tk  •  reaaining  cross  sections  ere  aost  easily  eoapnted  by 
integrating  the  norael  ooaponent  of  the  Poyatiai  vector  over  the 
fiber  surface,  neglecting  end  surface  contributions.  For  the 
scattered  power  one  has 


I,  -  |  fte  f  dS  •  (E  x  B*).  (1-14) 

•  2  J  ~S  ~S 

Ve  are  interested  only  in  effects  due  to  the  x-ooapoaent  of  E,, 
which  is  assuaed  independent  of  axiauth  t  because  ka  <<  1. 
Consequently,  only  the  f-coaponeat  of  Bs  is  needed,  and  because 
the  inoident  aagnetic  field  H*  varies  as  cos  #  it  will  not 

see 

oontribute  to  the  aurfaoe  integral.  We  oan  thus  replace  B,  by  H 

mm 

"  Hi  +  Hs.  If  Eg  is  also  replaced  by  E  -  Ej,  then  Eq.  (1-14) 

«V  mm  a* 

beooae  s 

I,  -  |  le  1  dS  •  (E  x  H*>  -  \  Ee  I  dS  •  (B*  x  H*) .  (1-15) 

The  first  tern  on  the  right  baud  side  gives  total  power  flow  into 
the  fiber,  and  hence  is  the  negative  of  the  absorbed  power.  The 
second  tern  is  proportional  to  the  real  part  of  the  forward 
aaplitude,  or  the  extinction  croas  section  Dividing  each 

tera  by  the  incident  energy  flux  density  Bq  /2q0  to  get  cross 
sections,  Bq.  (1-15)  beooaea 


oe  -  (4n/k2)  Re  S  (9*,  6*)  -  o,  +  c.,  (1-16) 


vkiik  it  tk*  optical  tkiom  i4«ttifyU|  tki  extinction  oton 
section  with  the  ita  of  the  scattering  tad  absorption  «tou 
MttiOll. 


Proa  Bq.  (1-14),  again  ambatitnting  fox  I.  ns  bofoxo  (bat 

*m 

Ul  »*)  ia<  normalising,  one  gets 


Mow 


dS 


*  BaJ*d8  •  (I|  x  I*). 


(Bg  x  H*)  -  (a  dtdx )  (Bg)  Xp* 

mm  *0  • 


wkara  (Eg) a  is  given  by  the  integral  term  of  Bq.  (1-1),  and 
Hp*  -  (l/2xa)  I*(x) . 

Patting  everything  together,  and  making  aie  of  Bq.  (1-t), 
the  result  is: 


*,  • 


-4n  la 
k1 


a 


One  proeeede  similarly  with  the  first  term  on  the  right  hand 
side  of  Eq.  (1-13)  to  obtain  the  absorption  cross  section,  noting 
that  Bs  -  ijg  I(x)  .  to  get 


Finally,  in  order  to  show  explicitly  that  the  extinction 


-17) 


(1-18) 


cross 


section  is  the  sum  of  these  last  two  expressions,  note  that  the 
integrand  in  the  last  tera  of  Bq.  (1-13)  aay  be  written 


d»  •  (*l  a  ■*)  -  («  d#ds)  Ea  aim  eA  oikl  *°*  6i  (1/ 2ftft) I*  (*). 

IV  W  M 

F«m  l|.  (1-9).  tko  *- integration  la  im  to  itvolvi  only  tko  *0,  sad 
vitk  minor  effort  o»«  giti 


wklok  ky  iaofootlom  is  tko  su  of  Eqs.  (1-17)  ftftd  (1-1*). 

Fee  loag  fibere,  a  mftofal  okook  in  provided  ky  «0ft|iUa|  tko 

•baoxytloa  oron  itotloa  resetting  (too  tko  laflalto  flkor 

•olatloa.  Afftla  latogratlag  tko  total  field  Poyatiag  veoior  over 

tko  lateral  ftmxfaeo  4xak  of  tko  flkor,  aslag  asset  field 
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fmaatltloe  a«  glvoa  by  Wait,  tko  rotfttioaally  symmetric  portioa 
of  tko  field  eoatribatos 

•*  ■  ”  *o  I  A0/Ei  I2  to  l(in/«|i0)  (k»)  Ji  («•)]  4a»k  (1-10*) 

vkoro  la  tko  limit  of  iaterest  to  as 

Ba  sia  #A/ A0  —  J0  (a*) 

ka  <<  1 

♦  sia*  0A  [la  ( ka  sia  0i/2)  +  y  +  i«/2]  asJi  (as)  (l-20b) 

vltk  y  ■  0.3772  ...  (Baler's  ooastsat).  Comparisons  asiag  tkoso 


formal**  will  bo  skova  bolov 


i 


At  low  intuxlM  for  whisk  fibsr  ilsisiiox  are  mill 
teafttii  wit k  viTiltKjtk.  11  ilickieititii  triitmt  is 
it)»»rtiiti.  Tks  elsotr.o  fisld  B  Uiiii  tks  fibsr,  regarded  as 
a  Bayleigk  ftrtiols,  ai;  b»  written 

B  -  Bt  -  L4*P,  (1-21) 

wksrs  L  is  tks  ispsliritiii  tsotor,  P  is  tks  pol aria  a ti ot,  and  B* 
tks  kOBOgsneoas  applied  fisld  along  tks  axis.  In  tks  oais  of  a 
long,  tkin  spkeroid,  wkiok  skonld  bs  a  good  appr oxiaati  on  to  tks 
fibsr  in  tks  prsssnt  oiroustasoss.  tks  dspolarixing  faotor  is 


L  -  (a/h)*  tin  <2k/a)  -ll .  (1-22) 


Tks  polarisation  is  defined  as 

?  -  U/4n)  (s-l>  B,  (1-23) 


wkers  s  is  new  tks  relative  dieieotric  constant  (nornslised  by 
dividing  by  s0).  From  Bqi.  (1-21)  1  nd  (1-23)  ons  has  for  the 

internal  f i el  1 


B  - 


E, 


1  +  L  ( s-1) 


(1-24) 
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and  if  *  it  tikia  to  b«  eoaplex,  l.t.  ,  •-•*«■  is",  then 


U  -  B 


r  1  +  L  (t*-!)  -  ILs"  -| 
i  LTI  +  L  («•-!)  1*  +  (L»,,Tr“  J 


(1-23) 


Not*  that  ia  the  paasistatic  approxiaa  tioa  tkitt  it  botk  ta  in- 
pkit*  tad  •  qudntari  ooapoaent  of  tko  internal  field. 


Tk •  oloettio  dipole  aoaeat  of  tke  pertiole  it  defiaed  n 


P  -  PV-. 


(l-2«) 


where  Vp  it  partiolo  volaae;  tke  po . ar ia abi 1 i ty  a  it  then  defiaed 
in  terat  of  tke  dipole  aoaeat  by 


P  *  «  Hi  . 


(1-27) 


Proa  tke  preceding  eqtatioai  tke  deaired  ezpretsioa  for 
pol aria abi 1 ity  it  given  by 


4w 


(«*-l)  [1+L  (t'-lM  »  Ls"  2  ♦  is" 
[1+L  (e  '-1^  P  +  (La"  T* 


( 1-2  8) 


Tke  absorption  and  scattering  cross  sections  are  siaply  expressed 
in  terns  of  a;  one  has,  respectively. 


o a  -  4xk  Ia  a. 


a*  -  (8n/3)  k*  lal2 


(1-29) 


(1-30) 


I 
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Thetc  eqvetioaa  kiv«  good  theoretical  jaatif ieation  within 
thi  tango  kh  &  a/2.  \Vhea  kk  ozeoodo  this  linit.  tko  effective 
wain#  of  tko  depolarizing  faotor  in  Bq.  (1-22)  it  redaoed.  This 
la  kooaaao  tko  dipolo  field  dao  to  tko  iadaoad  ekargot  at  tko 
ondo  of  tko  partielo  ia  no  longot  etrictly  a  ont  of  pkaao  with 
tko  applied  field.  In  order  to  take  tkio  diacrepaaoy  iato 
aoooaat,  wo  have  deviaed  a  nodifiod  dtpolarixiag  faotor  given  by 


L  « 

\ih  '  +  (H)  ^ 

]  fu  (f.  -1  1 

J  L .  *  J 

At  waveleagtha  for  wkiok  kk  >  a,  tko  depolarizing  faotor 
iaereaeet  at  k* .  Tkit  ia  oonaitteat  witk  oar  roaaoning  regarding 
tko  kigk  froqaoaoy  behavior  of  tko  iaookoronoo  of  tko 
depolarizing  field.  It  ia  beeaate  of  tko  inolmtion  of  Bq.  (1-11) 
tkat  wo  refer  to  tkit  tkeory  aa  oar  Bztondod  Qaaeietatio  (BOS) 
tkoory.  Vo  tkoald  note  tkat  tko  faotor  (1/a)  appoariag  in  Bq. 
(1~31)  vat  arrived  at  tkroagk  nany  oonparioont  between  the  BQS 
and  tko  Variational  tookniqao  diaeaotod  in  tko  provionv 
aabooo  tion. 

Altkoagk  tko  above  (BOS)  tkoory  ia  only  approxinate,  it  hat 
tko  very  atofal  proportioa  of  being  tinplo  and  anal  vtl  o.  It  oan, 
therefore,  bo  aood  to  provide  roaoonablo  approxinato  prediction! 
for  oatot  of  iatoroot.  Bowovor,  taok  prediction*  (oapccially  in 
tko  rogino  kk  >  1)  tkoald  bo  altinatoly  refined  naing  tko 
variational  tookniqao.  Cooperative  reaolts  will  be  shown 
aabaeqaently  in  tkit  report. 
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1.4  Tko 


1 


1b  order  to  iffljr  tk«  above  result*  to  the  iufrired  and 
visible  regions,  as  veil  as  tke  aicrowsve  case,  it  is  aeoessary 
to  build  is  tke  optical  behavior  of  tke  f  lbeji.  I*  this  aid  tbe 
folios  lag  seotioo  we  deal  with  (1)  tke  imelusioa  of  tk*  optical 
properties,  aid  (2)  imolasiom  of  tke  depemdeaoe  of  eleotrioal 
conductivity  upoa  partiele  radius  a  ad  eleotroa  aeaa  .free  patk.  ia 
botk  tke  Oalerkia  aad  quasistatie  aodels. 

la  a  reoeat  paper  by  Ordal  et  al.1^,  tke  applioatioa  of  tke 
Drade  aodel  for  tke  predietioa  of  toaplei  optical  dieleotrie 
eoaataat  was  compared  with  reasirad  values  of  tk*  real  aad 
imaginary  parts  of  tk*  optioal  dielectric  ooastaat  for  a  anaber 
of  aetals  (Al,  Cu,  Am.  Pb,  Ag.  aad  V).  Tabular  experimental 
results  vers  also  given  for  Pe.  Pt*  Co.  Ni,  Ti  aad  Pd.  Tkis 
aodel.  whiek  is  based  oa  the  free  eleotroa  theory  of  aetals,  is 
ia  surpri siagly  good  agreement  wltk  tk*  observed  experimental 
results.  Ve  realise  tkat,  for  certain  transition  eleacats  suoh 
as  Pe,  tk*  aodel  kas  drawbacks.  Por  euok  cases,  one  suit  resort 
to  tk*  us*  of  tabular  experiaeatal  data. 

Por  our  present  disoussioa,  we  okoos*  Cu  as  tk*  substance 
coaprisiag  our  fibers,  aad  will  utilise  tk*  Drud*  aodel  ia  tk* 
calculation  of  tk*  various  electromagnetic  cross  sections.  A 
good  exposition  of  tkis  is  given  in  Wooten'*  book,  in  wkich  the 
real  aad  imaginary  parts  of  the  dielectric  constant  take  tk*  fora 


IS 


•  '  -  1 


(1-32) 


l  + 


«  "  m 


*»  tl  ♦  (•*)*] 


(1-33) 


D»i*f  aki  salts,  tit  plasaa  fraquoaoy  wp  it  giraa  by 

•$  ■  naJ/i«0  (1-34) 


la  vkiok  a  -  alaotroa  daaaity  (a-3),  a  -  alaotroaio  changa,  a  - 
affaetiva  aaaa  of  tka  alaotroa,  aad  t0  -  paraittivity  of  fraa 
a pa  oa  *  ( 1/ 36  a)  a  10~9  farads/a. 


Tka  qaaatlty  x  la  tka  alaotroa  ralaxatioa  tiaa,  vkiok  it  .tka 
tiaa  raqalrad  for  raad  oa  is  a  t  loa  of  tka  aoaaataa  vaetor  of  an 
alaotroa  la  tka  (aatallio)  lattica.  For  oar  parpoaai,  it  la 
iaatruotlva  to  oast  tka  dialaotrio  constant  la  taraa  of  tka  lov 
freqnaaoy  alaotrioal  conductivity  a,  givaa  by 

«  -  aa2r/a  .  (1-35) 

Nov  ■  o/x i0,  aad  Bqa.  (1-32)  aad  (1-33)  aay  ba  vrittan  aa 


a  '  -  1 


£1  (  1 

•  o  V  1  +  (at)1 


(1-36) 


a 


•  •  - 


a 

•  a0  [1  +  (v»x )  ^  J 


(1-37) 
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Typical  values  of  the  relaxation  tine  t  are  on  the  order  of  10_14 


sec. 


Ik  is  easily  shorn  that  Eq .  (1-36)  can  be  written  in  terms 

of  b  "  : 


e  1  =  1  -  e  "  (at) 


(1-38) 


This  equation  shows  that,  for  all  frequencies  significantly  below 
the  visible  and  infrared  (i.e.,  the  microwave  region),  le'l  << 

I  8  "I  .  Also,  from  Eq.  (1-37),  we  see  that  a1'  goes  to  its  low 

frequency  value  6"=  a/ wa0  for  (<er)^  <<  1. 

The  reason  for  the  above  analysis  is  to  determine  whether  or 
not  the  Erode  model  can  be  utilized  at  1  ow  frequencies.  Although 

the  low  frequency  value  of  «'  differs  significantly  from  a  value 

of  unity,  which  is  normally  assumed  for  metals,  the  ratio 
I  e  "/ e  '  I  will  always  be  very  large  when  (ut)2  <<  i,  Therefore, 

the  use  of  the  Drude  model  throughout  the  region  10“^  m  <_  1 

10  1  m  appears  to  be  justified,  and  we  feel  confident  in  using 

Eq  s  .  (36)  and  (37)  in  the  derivation  of  the  electromagnetic  cress 
sections  throughout  this  entire  wavelength  range  for  appropriate 


ma  teri al s . 


1.5  Til  kedaesd  Conductivity 


W hen  one  or  aore  dimensions  of  a  conductive  material  (metal 
or  semiconductor)  are  on  the  order  of  the  mean  free  path  of  the 
conduction  electrons,  electron  collisions  vith  the  surface  will 
significantly  reduce  the  mean  free  time,  and  hence  the  mean  free 
path  A  given  by 

A  =  vpx  ,  ( 1~3  9 ) 

where  vp  =  Fermi  velocity  and  x  =  relaxation  time  discussed 
previously.  Since  the  electrical  conductivity  is  proportional  to 
x  (see  Eq.  -35),  the  conductivity  will  also  be  reduced. 

The  classic  work  on  this  subject  was  done  in  1938  by 

Fuchs.19  In  a  more  recent  paper,  Dingle  reviews  the  subject  and 

2  0 

provides  some  useful  nuaexiosl  computations.^  The  key  equation 
in  Dingle's  paper  is  his  Eq.  (2.3)  which  relates  the  effective 
conductivity  a  to  the  bulk  conductivity  a  Q,  as  a  function  of  mean 
free  path,  wire  radius  a  and  the  quantity  e  that  is  the 
probability  of  an  elastic  collision  at  the  surface: 

o  =  o0  U  -  |  (1  -  e)  (£)]  .  (1-40) 

A  value  of  e  =  1/2  is  frequently  used  as  the  surface  scat¬ 
tering  coefficient.  Taking  this,  we  have 

c  =  oc  [1  -  r|  <£)]  .  (1-41) 


We  have  had  difficulty  in  obtaining  numerical  values  for  *ean 
free  paths  f  roe  the  literature.  However,  Kittel21  provides  a 
good  background  as  well  as  quantitative  data  for  a  number  of 
■  etala.  Taking  copper  as  the  subject  material,  a  value  of  * 
4.2  x  10  ®m  is  given  in  Table  10.1  of  Kittel's  book.  Using  Eq. 


(1-39)  . 

and  taking  the 

Fermi 

velocity  vp  =  1.6  x  106 

m/  sec 

21 

one 

obtains 

a  mean  free  time  of 

—  14 

t  =  2.6x10  sec. 

This 

is 

in 

fairly 

good  agreement 

with 

the  accepted  value  x  ■* 

1.9  x 

10 

-14 

sec . 


Utilizing  the  above  value,  one  gets  for  copper 


r.  7.9  i  10  , 

a  =  oQ  [1  - - - - ] 


(1-42) 


where  the  radius  a  has  the  unit  of  meters.  This  equation  shows 

-  O 

that,  if  a  **  1.6  x  10  m  (160  Angstroms),  the  conductivity  is 
roughly  half  its  bulk  value.  If  the  radius  is  0.1  mioron,  the 
conductivity  is  92%  of  the  bulk  value.  This  exercise  was  done  to 
show  that,  indeed,  one  must  consider  the  particle  size  effect 
upon  electrical  conductivity,  when  competing  absorptive  and 
extinction  properties  of  thin  metallic  fibers  or  films. 


1.6  litiairitil  Kiitlti  fci  Sstfkit*  ta4  Itoa 


The  precede  &|  thioiy  hat  been  reduced  to  computer  codes 
including  grspbiosl  outputs,  and  some  numerical  results  are 
presented  here  for  tbe  two  cases  of  graphite  and  iron.  It  should 
be  noted  that,  as  part  of  our  AFOSE  contract,  we  have  received  a 
Hewlett  Packard  HP  90  00  Mod.  S20  oomputer.  which  is  employed  for 
all  computations. 

Considering  first  graphite,  the  sequence  of  operations  is  as 
follows:  Taking  a  bulk  conductivity  of  10$  mho/m  and  mean  free 

path  of  1.4  z  10"®  m,  the  reduced  conductivity  is  obtained  from 
Eq.  (1-41)  upon  specifying  the  fiber  radius.  Using  this  result 
and  further  assuming  an  electron  relaxation  time  of  1.4  z  10  15 
sec.,  the  complex  optical  dielectric  oonstant  is  obtained  as  a 
function  of  frequency  from  Eqs.  (1-36)  and  (1-37).  The  surface 
impedance  of  the  fiber  can  now  be  computed  from  Eq.  (1-2),  at 
which  point  the  scattering  computations  can  be  oarried  out. 

Some  typical  B-plane  differential  scattering  patterns  for 
graphite  are  shown  in  Figs.  1—2  thru  1-4.  The  fiber  has  radius 
of  one  micron  and  length  of  110  microns.  The  incident:  wavelength 
is  30  microns,  so  that  kh  *  12  and  the  fiber  is  several 

wavelengths  long.  The  angle  of  incidence  ranges  progressively 
from  broadside  to  near  end-on  aspect  in  the  three  figures,  and  as 
one  would  expect  a  fairly  large  main  lobe  is  seen  in  the  forward 
direction  in  all  cases,  with  much  smaller  side  lobes  in  other 
dire  ctions. 


20 


Figure  1-2.  Differential  scattering  pattern  ia  shown  va. 

scattering  angle  for  a  graphite  fiber  (noraal 
inside  nee ) . 


DIFFERENT I RL  SCATTERING 
Cross  Met  t  on/)  vs  Scatter  tng  angle 


e  36  66  96  120  150  160 


Figure  1-3.  Differential  scattering  pattern  is  shown  vs. 

scattering  angle  for  a  graphite  fiber  (incidence 
at  60  deg.  froa  the  axis). 


The  integrated  cron  eeetione  (ecattering,  absorption  and 
extinction)  are  shown  for  the  aane  graphite  fiber  in  Pig.  1-J  < 
for  wavelengths  ranging  fron  one  aioron  np  to  10  on.  It  should 
be  noted  that  we  have  done  a  planar**  averaging  (fiber  axis  in  the 
plane  forned  by  the  incident  and  observation  directions).  The 
fnll  random  orientation  results  are  obtained  by  redncing  the 
cmrves  in  the  figures  by  an  additional  3  db.  Note  that  the 
extinction  cross  section  is  precisely  equal  to  the  sun  of  the 
eoattering  and  absorption  cross  sections  at  all  wavelengths,  as 
discussed  earlier.  Low  frequency  Rayleigh  behavior,  (l/X.)*  for 
the  scattering  and  (1/X)2  for  the  absorption,  is  clearly  evident 
at  the  longer  wavelengths.  At  short  wavelengths,  the  scattering 
cross  section  becomes  dominant. 

The  corresponding  integrated  cross  sections  for  a  graphite 
fiber  slightly  larger  in  diameter  and  about  ten  times  longer  are 
shown  in  Fig.  1-6  (see  figure  for  precise  dimensions).  As  one 
would  expeot,  the  peak  values  of  all  the  cross  sactions  now  occur 
at  wavelengths  ten  times  greater,  although  the  magnitudes  of  the 
peaks  are  virtually  unchanged. 

The  different  methods  of  computing  absorption  cross  section 
are  oompared  in  Fig.  1-7  for  a  small  graphite  fiber  (radius  * 
0.02  microns,  length  =  16  nicrons).  The  Galerkin  (or  varia¬ 

tional)  results  correspond  to  the  curve  with  oscillations  at  the 
shorter  wavelength  end  of  the  spectrum,  similar  to  the  results 
seen  earlier.  The  extended  quasistatic  approximation  is  seen  to 
match  the  variational  curve  identically  at  the  lower  frequencies 
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logfCrfii  Sset/'untt  *tl) 


up  to  about  kt  ■  2x  (a,  -  1.6  x  1 0“ -  ■)  tad  oontlnnti  to  folio* 
tilt*  voll  at  higher  (r iqwioiti.  Tba  high  frequency  infinite 
cylinder  appr oxiuati oa  (ihovs  beginning  at  k£  *  1)  also  gives  a 
good  fit  to  tba  variational  osrva,  al though  it  it  slightly 
higher. 

Tha  oorr aspoadiag  curves  ara  shown  in  Fig.  1-8  for  tha  sane 
fiber  whose  length  has  bean  inoraasad  by  a  faotor  of  eight. 
Kesults  ara  coupletely  analogous  to  those  of  Fig.  1-7,  although 
of  course  tha  layleigh  region  has  no*  aoved  offsoala  to  tha 
right,  and  high  frequency  oscillations  of  tha  variational  curve 
ara  nueh  sore  rapid. 

Turning  no*  to  tha  oasa  of  iron  fibers,  one  has  bulk 
conductivity  of  10?  ihot/a  and  naan  free  path  of  3  x  10-11  m. 
From  Bq.  (1-41)  one  seas  that  tha  bulk  conductivity  requires  no 
correction  this  tine,  even  for  wire  radii  as  snail  as  10"*  n. 
Choosing  an  alaotron  relaxation  tine  of  4.0  x  10”13  *#©. ,  the 

conplex  dielectric  constant  is  again  obtained  as  a  function  of 
frequency  f  ron  Eqs.  (1-36)  and  (1-37). 

Tha  differential  scattering  patterns  for  iron  are  shown  in 
Figs.  1-9  thru  1-11  for  broadside,  (0  and  30  deg.  incidence, 
respectively.  Conparing  with  the  earlier  Figs.  1-2  thru  1-4  for 
graphite,  the  angular  distributions  are  seen  to  be  quite  sinilar. 
The  peak  cross  sections  per  square  wavelength  are  seen  to  be 
snaller  fox  iron  by  about  2  1/2  orders  of  nagnitude  (conpare  the 
respective  values  of  Y  in  the  figures),  but  on  the  other  hand. 
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the  aotual  f ib«r  woluae  for  iron  hot  been  reduead  by  about  5  1/2 
ordora  of  magnitude,  ao  that  on  a  par  unit  voUai  baaii  Iroa 
aoattara  amok  aora  of  f  ioiently,  aa  oaa  would  axpaot. 

Oriaata tioa-awarafad  eroaa  aaotioaa  for  Iroa  fibara  ara 
akowa  ia  Figa.  1-12  tkru  1-14,  for  progr eaaiwely  loagar  fibara 
agaia  witk  radiua  .01  aieroaa.  For  tka  akortaat  fibar,  laagtk  “ 
S  aieroaa*  akowa  ia  Fig.  1-12,  tka  abaorptioa  oroaa  aactioa  ia 
aooa  to  riaa  wa.  frequency  aad  paak  fairly  akarply  ia  tka 
wioiaity  of  k t  »  n.  Tka  aoattariag  oroaa  aaotioa  ia  quite 
aagligibla  ia  ooapariaoa,  axoapt  at  tka  akortaat  wawelengtha,  aad 
ooaaaquaatly  tka  abaorptioa  aad  aztiaotioa  omrwea  ara  praotioally 
iadi a tiagui akabl a .  Tka  daakad  ourwaa  giwa  tka  quaaiatatio 
appr oxiaati oaa,  wkiok  ara  aaaa  to  ba  excellent  for  k L  ^  a.  Tka 
apptoxlaatc  abaorptioa  oroaa  aaotioa  aatokaa  tka  wariatioaal 
ooiaputatioa  wary  wall  up  to  about  k l  -  n,  aad  ia  aoaowkat  kigkar 
at  tka  akorter  wawalaagtka.  Siailar  ooaaaata  apply  for  tka 
approxiaata  aoattariag  oroaa  aaotioa,  al  tkougk  tka  lattar  kaa 
baaa  out  off  aad  ia  aot  axpaotad  to  apply  bayoad  k£  »  x. 

Reaulta  ara  aiailar  wkea  tka  fibar  laagtk  ia  doubled,  aa 
akowa  ia  Fig.  1-13,  although  the  paakiag  of  the  abaorptioa  (and 
aztiaotioa)  oroaa  aaotioa  ia  aot  quite  ao  proaouaead.  Vhaa  the 
fibar  length  ia  iaereaaed  by  a  further  order  of  uagaitudo  (Fig. 
1-14),  the  riaa  to  paak  falsa  a  ia  aaaa  to  occur  at  auch  longer 

■  aval  aaafk*  1  AKA  ■  A  HI  1  1  A  f  n  f  T1  W  ft  a  Afltifiiflkti.  ft  fid  t  H  fl 


abaorptioa  (and  a xtiac tioa )  oroaa  auction  now  showa  a  broad 
plateau  extending  over  aora  than  a  decade  in  frequency. 
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It  it  i  •  t  •  r  t  i  t  i  a|  that,  nidi  (tot  the  ton  rapid 
oioillatioti  occurring  for  thi  longer  fibart  ia  the  abort 
vaveleugth  region,  due  to  partial  dependence  oa  trigo  oaatric 
faaotioaa  of  argaaaat  k£,  tba  oroaa  Motion  ara  effectively 
iadapaadaat  of  fibar  laagtb  at  tba  abortaat  wavelengths. 

Soai  aatpla  eoapvtati oni  of  abaorptioa  oroaa  aaotioa  for 
iroa  ara  shovn  ia  Figa.  1-15  aad  1-16.  Tba  fibar  radii  ara  2  and 
8  aioroat,  raipa  o tivaly ,  aad  tba  laagtb  ia  64  aiorona  ia  both 
oaaaa.  Tba  abaorptioa  ooapttad  by  tba  variational  aatbod  ia 
aiailar  to  r'  salti  abova,  shoving  a  paak  at  a  wavelength  of  about 
100  aiiorona.  In  Fig.  1-15,  tba  axtandad  qnaaiatatio  ourva  tracks 
tba  variational  raailta  fairly  vail  azoept  at  tba  shorter 
vavalangtba,  vbara  it  dees  not  daereaae  aa  it  should  (tba  raaaon 
for  this  diaorapanoy  ia  not  praaantly  understood).  Tba  infinite 
oylindar  approxiaa tion,  on  tba  otbar  band,  doaa  raaaonably  wall 
at  tba  sbortar  vavalangtba.  For  the  larger  fibar  radios  of  Fig. 
1-16,  both  of  tba  approsiaata  computations  ara  aaan  to  be  in 
raaaonably  good  agraaaant  vitb  tba  variational  result. 


It  ia  straight! orvard  to  ooaputa  tba  cloud  aass  par  square 
aatar  required  to  produce  a  20  db  target  signature  reduction, 
using  either  graphite  or  iroa  particles.  The  signal  decays 
exponentially  vithin  the  cloud,  and  a  10  db  reduction  on  each  lag 
of  tba  round-trip  path  requires  that 


-  n3  l 

ext  «  _  . 

a  •  0 . 1 ,  or 


n8extLc  "  2’3* 


(1-43) 
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where  n  is  the  numbe  r  of  fibers  per  uni  t  volume  end  Lc  is  path 
length  through  the  cloud.  We  include  an  additional  fact'  >-  of  1/2 
in  the  extinction  cross  sections  9  given  earlier,  so  that  all 
polarization  and  (random)  orientation  affects  are  now  completely 
accounted  for.  Strictly  speaking,  Eq.  (1-43)  aasur.es  that 
absorption  dominates  over  scattering  for  the  single  particle.  If 
this  is  not  the  case  then  multiple  scattering  effects  must  be 
core  carefully  accounted  for,  using  the  radiative  transfer 
analysis  to  be  discussed  shortly. 

For  a  section  of  the  cloud,  having  cross-sectional  area  Ac 
and  length  Lc,  the  total  mass  of  fibers  N  is  given  by 

M  -  AcLc  hVF  pp,  (1-44) 

where  Vp  is  the  volume  of  a  single  fiber  and  pp  its  density.  For 
the  case  at  hand,  setting  Ac  ®  1  and  substituting  in  from  Eq . 

(1-43)  gives 

M  -  2.3  Pp/Oext/Vp)  .  (1-45) 

Using  this  result,  from  the  extinction  cross  section  por  unit 
volume  data  of  Fig.  1-5,  and  taking  a  density  4  gm/cm^  for 
graphite,  one  obtains  the  mass  requirod  at  various  wavelengths  as 
listed  in  Table  1-1.  Note  that  the  values  in  parentheses  may  be 
underestimates  because  the  fibers  are  not  primarily  absorbers  at 
those  wavelengths  (see  Fig.  1-5). 
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8  i  gni  £  icantly  let*  am  is  required  eaploying  iron  fibers, 
es  can  be  seen  froe  the  last  solses  of  the  table.  Sere  we 
employed  the  data  of  Fig.  1-14,  assuaing  -  7  ga/ ca3  for  iron. 
Note  that  the  increase  in  density  is  aore  than  ooapensated  by  the 
higher  conductivity ,  which  in  turn  allows  one  to  eaploy  a  greatly 
reduced  fiber  radius. 


Table  1-1.  The  aass  of  graphite  or  iron  fibers  required  per 

square  aster  of  cloud  cross  section  for  a  20  db 
reduction  in  signal  strength  (values  in 
parenthesis  aay  be  undere  a  tiaate  s  ;  see  text). 


enath  k  ( a) 

Hass  of  arauhite  (an) 

Hass  of  iron  (aa) 

1 

o 

H 

(180) 

(1.8) 

•n 

1 

O 

(35) 

0  .079 

io-4 

(6  .3) 

0.027 

io-3 

72 

0.025 

M 

1 

O 

H 

7500 

0  .16 

10_i 

— 

10 
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21. 7  UlittiT*  Tramsfer  ia  Pib«r  Clouds 


In  order  to  ooapnte  the  tren mission  end  reflection  by  e 
cloud  of  fibers,  including  multiple  scettering  effects,  we  enploy 

radiative  transfer  techniques,  using  e  ooaputer  progrea  developed 

2  2  2  3 

earlier.  *  Essentially,  the  computation  employs  the  van  de 

Hulst  doubling  method  in  a  slab  geometry,  with  Taylor  series 
r  e  p  r  e  a  e  n  t  a  t  i  o-n  for  the  thin-layer  starting  values  of  the 
transmission  and  reflection  matrices. 

For  a  cloud  of  identical  fibers  randomly  oriented,  as 
considered  here,  the  required  inputs  from  the  single-scattering 
results  are  the  albedo  (ratio  of  scattering  to  extinction  cross 
section)  and  the  differential  scattering  pattern  of  the  single 
fiber.  The  calculation  has  been  carried  out  for  both  graphite 
and  iron,  assuming  monochromatic  radiation  normally  incident  on  a 
slab  region  containing  the  fibers.  In  each  case  both  wavelength 
and  fiber  dimensions  are  chosen  so  that  the  scattering  and 
absorption  cross  sections  are  roughly  ocmparable.  Results  will 
be  given  fox  seleoted  values  of  the  optical  depth  of  the  slab 
(optical  depth  is  obtained  by  multiplying  the  thiokness  by  the 
number  of  fibers  per  unit  volume  and  the  extinction  cross  section 
of  one) . 

Consider  first  graphite.  For  an  incident  wavelength  of  30 
microns  and  particle  dimensions  as  listed  in  Fig.  1-17,  the 
fractions  of  incident  energy  reflected  from  the  front 
(illuminated)  face  and  transmitted  through  the  bank  (shadowed) 


fie*  m  shows  in  thi  figure  at  a  function  of  |i  ■  001  6,  where  0 
it  tkt  loitttriti  angle  measured  frot  thi  original  incidtnt 
direction.  Optioal  depth  it  0.5  in  this  iaataace.  Note  firat 
the  4eltt-f emotion  behavior  (trnnoated  at  unity  for  clarity)  of 
the  traasaitted  power  in  the  forward  direction  |i  **  1 .  Thia 
repretanta  the  coherently  trananitted  tern  which  hat  been  neither 
acattered  or  abaorbed.  Aaide  iron  tone  ninor  f  luctua  ti  one,  the 
balaaoe  of  the  trananitted  flux  it  teen  to  be  eaaentially 
isotropic,  at  it  the  reflected  flux,  although  both  tend  to  be 
alightly  higher  in  the  graxing  direction*  ( |i  »  0)  . 

Figure*  1-18  thru  1-20  give  the  corresponding  results  with 
sueceaive  doubling*  of  the  slab  thickness,  holding  all  other 
parameter*  fixed.  The  trananitted  flux  is  teen  to  gradually 
deereaae,  at  one  would  expect*  and  in  addition,  bee  one 
ooneentrated  nor*  in  the  forward  direction.  Note  that  the  delta- 
function  contribution  ia  difficult  to  ascertain  graphically. 
Thia  is  no  problen*  however*  beoaus*  we  know  apriori  that  the 
fraction  of  coherently  trananitted  flux  is  given  sinply  by  e~T 
(t  ”  optioal  depth).  The  reflected  flux  distribution  is  seen  to 
beoone  constant  at  the  largest  optioal  depths;  thi*  it  of  oourae* 
just  the  distribution  one  would  obtain  f row  an  infinitely  thick 
slab. 


These  curves  have  been  integrated  to  obtain  the  total 
fractions  of  incident  flux  reflected  and  trananitted,  and  the 
result*  are  listed  in  Table  1-2.  In  each  case*  of  coarse,  the 
balance  of  the  incident  flux  i*  absorbed  within  the  fibers.  The 
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Figure  1-17 


The  refleoted  end  transmitted  intensities  ere 
plotted  ts.  oos  0  (6  is  the  soattexing  angle)  for 
a  olond  of  graphite  fibers  of  optioal  depth  0.5. 
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Figure  1-18. 


The  reflected  and  transaitted  intensities  are 
plotted  vs.  cos  8  (8  is  the  scattering  angle)  for 
a  cloud  of  graphite  fibers  of  optical  depth  1. 


Tk  •  triotiom  of  (integrated)  reflected  and 
traaeaitted  flux  act  given  vi.  optical  depth  for  a 
cloud  of  graphite  fibara  (albedo  ■  0.65). 


. . 
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Table  1-2 . 


ftjfcliJUQ  ff.tatfc 

tef lcctad  flux 

iill 

0.5 

0  .083 

0.729 

1  .0 

0.120 

0.521 

2  .0 

0  .145 

0  .257 

4.0 

0.151 

0.058 

8.0 


0  .152 


0.003 


total  reflected  portion  it  aeen  to  aatnrate  at  about  15%  for  the 
thicker  alaba.  It  ia  al  ao  intereating  to  note  that  for  the 
thiokeat  alab  the  coherent  tern  nakea  only  a  nepigible 
contribution  to  the  total  tranaaitted  flux. 


The  oorreaponding  reanlta  for  a  cloud  of  iron  fibera  arc 
ahovn  in  Figa.  1-21  through  1-24,  thia  tine  uaing  an  inoident 
wavelength  of  1.1  aicrona  and  appropriately  rednoed  fiber 
diaenaiona,  aa  liated  in  the  figurea.  Moat  of  the  coaaenta  aade 
above  in  connection  with  graphite  are  aeen  to  apply  here,  alto, 
but  there  are  two  aignificant  change  a.  Firat,  the  reflected  flux 
pattern  ia  aeen  to  be  aoaewhat  weaker  and  haa  already  beooae 
effectively  independent  of  alab  thickneaa  at  an  optical  depth  of 
0.5.  Second,  the  traneaitted  flux  falle  off  nuch  faater  with 
inoreaaing  thiokneaa.  Both  of  theae  ohangea  axe  attributable  to 
the  anallex  albedo  encountered,  whioh  ia  of  courae,  juat  another 
way  of  aaying  that  abaorption  haa  beooae  dominant,  over  acattering 


in  thia  caae. 


CyllntrUli  Baaltarara 
Reflection  T  r  aneml ee  1  on 


Fl|«t«  1-21.  Thi  rtf  1 eotid  and  trananittad  Intanaitiaa  ara 
plotted  va.  ttoa  0  (8  ia  tha  aoattaring  angle)  for 
a  dead  of  iron  fibere  of  optioal  dapth  O.S. 


Cylindrical  Scattarara 
Reflection  Trenemi salon 


-1  -.5  0  .5  l 


Figure  1-22.  The  reflected  and  trananittad  intanaitiaa  ara 
plotted  vb.  coa  8  (6  ia  the  aoattaring  angla)  for 
a  cloud  of  iron  fibara  of  optical  dapth  1. 


-.5  0  .5 


reflected  md  tnsiaitted 
tted  v«.  cos  8  (6  is  ths  scst 
loud  of  iron  fibers  of  optiosl 
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Tkt  eottiiyo>4ia|  i>t«|nttd  {Uin  for  iroa  art  liitid  ▼«. 
oftioal  4*ftk  ia  Tabla  1-3.  Oaa  aaaa  fro*  tka  tabla  tkat  tka 
total  raflaotad  flax  iadaad  okaagaa  taty  littla  witk  lmoraaaiag 
alak  tkiokaaaa,  aad  at  tka  aaaa  tiaa  tka  traaaaittad  flax  ia 
aoaovkat  aaallar  a  ad  falliag  off  aora  rapidly  tkaa  ia  tka  oarliar 


Tabla  1-3. 


Tka  fraotioaa  of  (iatagratad)  raflaotad  and 
traaaaittad  flax  ara  giraa  ▼*.  optioal  daptk  for  a 
eload  of  iroa  fibart  (albado  ■  0.332). 
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SECTION  II .  THERMAL  RADIATION  BY  SMALL  PARTICLES 


2.1  Introduction 

In  tho  calculation  of  electromagnetic  scattering  and  absorption  by  snail 
particles  having  at  least  one  dimension  on  the  order  of  the  wavelength  of  the 
incident  radiation,  it  Is  not  unconaon  that  the  absorption  (as  well  as  the 
scattering)  cross  section  can  be  significantly  greater  than  the  projected  area 
of  such  particles.  This  is  easily  understood  from  the  electromagnetics  point 
of  view  and  causes  no  serious  concern.  However,  if  we  consider  the  thermal 
radiation  from  such  a  particle,  <me  interesting  questions  arise  with  respect 
to  its  thermal  emissivity,  its  thermal  radiation  spectrum,  and  the  connection 
between  such  radiation  and  classical  radiation  theory.  The  purpose  of  this 
paper  is  to  address  these  questions  and  to  propose  a  quantitative  model  for 
the  calculation  of  the  electromagnetic  radiation  spectra  of  small  absorbing 
particles.  Of  particular  interest  will  be  particles  of  the  type  discussed  in 
Section  I,  namely  thin  conductive  cylindrical  particles  having  large 
absorption  efficiency  factors. 


2.2  Classical  Radiation  Theory 


In  this  section  we  summarize  those  aspects  of  the  classical  theory  of 
radiation  that  are  pertinent  to  the  discussion  to  follow. 


Thermal  Equilibrium.  Take  the  classical  case  of  a  perfectly  insulated 
spherical  cavity  whose  temperature  is  T  and  whose  wall  absorbs  100%  of  all 
radiation  incident  upon  it.  If  all  elements  of  the  wall  are  to  be  in  thermal 
equilibrium,  thin  each  element  must  radiate  100%  of  the  radiation  incident 
upon  that  element.  Since  this  must  be  true  at  all  temperatures,  it  follows 
that  the  absorptivity  a  and  the  emissivity  e  of  the  wall  must  each  be  equal  to 
unity.  The  resulting  radiation  is  well  known  as  blackbody  radiation,  whose 
spectral  distribution  is  given  by  the  Planck  Radiation  Function. 

To  summarize, 


“bb  "  ebb  “  1' 


(2-1) 


where  ^  is  taken  to  represent  blackbody  properties. 


The  Planck  Radiation  Function  is  given  by: 


■*r 


Wbb(A,T) 


2wc2h 

A5(ehc/AkT  .  1}* 


(2-2) 


where  c-speed  of  light,  h-Planck's  constant,  k-Boltzmann ' s  constant, 
A-wavelength ,  and  T-absolute  temperature.  With  the  insertion  of  an 
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appropriate  constant  (10*^)  ,  the  units  of  Wbb  are  the  usual 
o 

watts/cm  -micron,  assuming  that  Eq.  (2-2)  is  originally  evaluated  in  MKS 
units . 

The  intensity  of  the  radiation  emitted  by  a  blackbody  is  the  integral  over  all 
wavelengths  of  the  Planck  Radiation  Function: 


Wbb(X,T)dX. 


(2-3) 


The  total  power  radiated  by  the  blackbody  is  obviously  the  product  of  the 
radiant  intensity  and  the  area  A^  of  the  blackbody: 


Pbb(rad)  “  Wbb  watts’  (2_4) 


Kirchoff's  Law.  Next,  place  within  the  cavity  a  body  having  area  A  and 
wavelength  dependent  absorptivity  a(A)  and  emissivity  «(A).  Since  the  radiant 
flux  is  uniform  within  the  cavity,  the  power  absorbed  by  the  body  is  given  by 


Pabs  “  A 


«(A)Wbb(A,T)dA. 


(2-5) 


which  are  that  these  quantities  are  the  ratios  of  absorptance  and  emittance, 
respectively,  to  those  of  a  blackbody. 
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Note  that,  in  Eq.  (2-5),  we  have  used  the  same  temperature  T  as  for  the 
(blackbody)  walls  of  the  cavity.  This  is  essentially  the  definition  of 
thermal  equilibrium,  namely  that  two  bodies  in  thermal  equilibrium  must  be  at 
the  same  temperature . 

The  power  radiated  by  the  body  is 


Prad  “  Aje(A)Wbb(AlT)dA-  (2‘6) 
o 

The  important  consequence  of  the  above  discussion  is  that,  if  the  body  is  in 
thermal  equilibrium  with  its  surroundings,  its  absorptivity  and  emissivity 
must  be  equal  at  all  wavelengths.  This  is  true  because,  in  thermal 
equilibrium,  Eq.'s  (2-5)  and  (2-6)  must  be  equal  for  any  value  of  temperature 
T. 


Therefore , 


a(A)  “  «<A). 


(7) 


Equation  (2-7)  is  what  we  might  term  the  "detailed"  statement  of  Kirchoff's 
law  of  thermal  radiation.  Ono  usually  sees  this  in  the  form  of  the  average 
values  of  a  and  e.  Since  (2-7)  is  crucial  to  the  subsequent  treatment,  it  was 
developed  in  what  we  believe  to  be  a  simple  but  fundamental  manner. 


2.3  Microscopic  Particles 


Absorbed  power.  The  power  absorbed  by  a  particle  at  any  specified  wavelength 

1  is,  by  the  definition  of  the  absorption  cross  section,  the  product  of  the 

incident  radiant  intensity  and  the  absorption  cross  section.  If  wavelength  of 

the  incident  radiation  is  X .  then 

o 


ab  s 


(2-8) 


In  the  more  general  case,  if  the  incident  radiation  occurs  over  a  spectrum  of 
wavelengths,  we  must  use  the  integral  form  of  (2-8): 


00 

-! 


Pabs  “  Kbs<A>Winc<A>dA’ 


(2-9) 


where  W^nc(A)  is  the  intensity  spectrum  of  the  incident  radiation. 

The  absorption  efficiency  factor  is  defined  as  the  ratio  of  the  absorption 
cross  section  of  a  body  to  its  projected  area  ApZ 

WA>  -  WA>/V  (2*10) 

Rewriting  Eq.  (2-9), 

00 

Pabs  -  ApJ<Wi>Wlnc<i>'U'  <2'U> 
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Note  that  Qa^s(^)  for  the  particle  in  Eq.  (2-11)  replaces  a(A)  for  the 
macroscopic  body  in  Eq.  (2-5).  We  write  explicitly, 

a(A)  “  Qabs(A)-  (2'12) 

Since,  in  Eq.  (2-7),  we  have  shown  that  a(A)  -  e(A)  for  any  body,  it 
directly  follows  that 


e(A)  “  Qabs(A)*  <2-13) 


Radiated  Power.  The  power  radiated  by  the  particle  is  therefore  given  by 

co 

fr.d  -  <2'14> 

o 

where  A-  total  area  of  the  particle. 

In  the  steady  state,  the  power  radiated  by  the  particle  is  equal  to  the  power 
absorbed.  Therefore,  Eq.'s  (2-11)  and  (2-14)  must  be  equal.  Equating  these 
provides  the  expression  for  total  radiated  power  by  the  particle,  as  a 
function  of  the  incident  radiation  spectrum: 


-4 


ao 


Prad  "  \  W^inc^**  "  Ap  WA>Wbb<'' 'T>dA  •  <2'15> 
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Using  the  expression  (2-2)  for  W^i  we  obtain  an  expression  involving  only 
the  temperature  of  the  particle  ,  its  geometry  (A  and  A^),  and  the 
integrated  spectrum  of  the  incident  radiation: 


00 


2irc2h 


00 


where  we  have  used  T@£f  In  place  of  T,  as  discussed  in  the  following 
section. 


Equation  (2-15)  has  the  units  of  intensity.  In  the  case  of  monochromatic 
incident  radiation  at  wavelength  XQ,  we  can  write: 


00 

4* 


iwc2h 


abs(^^/^rir  "  VW  V  WV  *  <2'17> 


Equations  (2-16)  and  (2-17)  can  be  used  to  compute  the  radiating  temperature 
of  the  body,  as  a  function  of  its  wavelength  dependent  absorption  efficiency 
factor  and  the  spectral  content  of  the  radiation  incident  upon  It.  Also,  the 
Integrand  of  these  equations  is  the  thermal  radiation  spectrum  of  the  body. 
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2.4  Interpretation 


Emissivity.  Equations  (2-16)  and  (2-17)  are  valid,  regardless  of  the 
numerical  values  Qa^a*  In  particular,  as  we  have  shown,  when  Qa^s  is 
larger  than  unity,  the  enissivity  is  correspondingly  large.  If  one  considers 
this  from  the  thermodynamics  point  of  view,  some  questions  arise,  such  as  "How 
can  a  radiating  body  have  an  enissivity  greater  than  that  of  a  blackbody?"  The 
following  is  a  brief  discussion  of  why  this  is  not  an  unreasonable 
circumstance . 

In  the  steady  state,  the  particles  comprising  the  body  are  elevated  to  higher 
energy  states  by  virtue  of  absorption  of  power  from  the  incident  wave.  The 
body  loses  energy  at  this  same  rate  by  virtue  of  the  (presumably)  random 
emission  of  photons  as  the  excited  particles  drop  into  lower  energy  states. 
In  the  case  of  a  'macroscopic'  body,  all  of  whose  dimensions  are  very  large 
compared  with  all  wavelengths  of  Importance,  it  is  easy  to  understand  how  the 
emissivity  must  be  dependent  only  upon  the  composition  of  the  body.  However, 
in  the  case  of  the  'microscopic'  body,  its  geometry  also  plays  an  important 
role  in  the  wavelength  dependent  probability  of  emission  of  photons.  Why 
should  this  be  so? 

The  answer  is  that  the  same  geometrical  properties  that  cause  the  body  to 
absorb  electromagnetic  energy  at  a  rate  proportional  to  its  absorption  cross 
section,  similarly  influence  thermal  radiation  by  the  quantum  mechanical 
system  of  partic]es.  A  detailed  analysis  of  this  problem  would  involve  the 


full  quantum  mechanical  treatment  of  the  system  of  particles,  in  which  the 
geometry  of  the  body  would  necessarily  be  Included. 

If  wa  think  about  the  absorption  from  the  point  of  view  of  Maxwell's 
equations,  but  consider  the  radiation  from  the  quantum  mechanical  point  of 
view,  the  picture  can  be  quite  confusing.  The  particles  under  consideration 
here  behave  in  some  respects  like  macroscopic  bodies  and  in  other  respects 
like  quantum  mechanical  systems  having  absorption  and  emission  spectra 
determined  by  the  modal  constraints  (which  depend  upon  the  composition  as  well 
as  the  geometry)  of  the  specific  particle. 

From  the  quantum  viewpoint,  if  the  dimensions  of  the  body  were  all  on  the 
order  of  an  optical  (radiating)  wavelength,  only  certain  wavelengths  of 
absorption  and  emission  would  be  permitted.  As,  e.g.,  one  dimension  is 
increased,  the  modal  structure  will  become  more  complex.  At  a  sufficiently 
large  dimension,  the  transition  probabilities  overlap  and  the  body  exhibits  a 
continuum  in  its  emission  and  absorption  spectrum. 

The  case  of  coherent  incident  radiation  and  incoherent  thermally  emitted 
radiation  would  not  pose  a  problem  if  the  analysis  were  carried  out  at  a 
sufficiently  fundamental  level. 


Temperature.  The  temperature  appearing  in  Eq.'s  (2-16)  and  (2-17)  can  also 
lead  to  some  questions  regarding  its  precise  meaning.  We  believe  that  there 
is  no  question  about  its  meaning  when  we  consider  the  thermal  radiation  from 
the  body.  But,  is  this  the  "thermodynamic"  temperature  that  one  could  in 


principle  measure  with  a  thermometer?  This  question  is  analogous  to  that 
cited  above  with  respect  to  the  emissivity.  The  body  under  consideration  has 
some  properties  that  are  most  easily  considered  from  the  macroscopic 
thermodynamic  point  of  view  and  other  properties  that  must  be  thought  of  from 
the  quantum  point  of  view.  The  concept  of  "temperature"  arises  from  the 
assumption  that  one  is  dealing  with  a  body  of  sufficiently  large  dimensions 
that  the  use  of  macroscopic  thermodynamic  quantities  is  appropriate .  To  avoid 
possible  contradiction,  we  will  use  the  term  "effective  radiating  temperature" 
Teff  in  the  discussions  to  follow,  and  as  in  Eq.'s  (2-16) , (2-17) . 
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2 . 5  Computational  Procedure 


In  this  section,  we  present  computational  results  obtained  through  the  use  of 
the  foregoing  theoretical  treatment.  Note  that,  in  the  use  of  Eq.'s  (2*16) 
and  (2*17  > ,  one  would  normally  be  provided  with  that  information  needed  to 
reduce  the  RHS  of  (2-16)  or  (2-17)  to  a  number.  The  mathematical  problem  is 
to  find  that  value  of  T^ff  in  the  integral  that  satisfies  the  equation. 

The  procedure  In  the  simpler  case  of  monochromatic  incident  radiation  is  as 
follows : 

(1)  Compute  Qabs(*0)  over  a  wavelength  range  sufficiently  broad  to 
include  \Q,  as  well  as  that  of  the  thermal  radiation  by  the 
particle . 

(2)  Evaluate  the  RHS  of  Eq.  (2-17). 

(3)  Perform  successive  numerical  integrations  of  the  LHS  of  (2-17)  to 
obtain  a  table  and  graph  of  Tfi£^  vs  I^nc(Ao) ■  The  wavelength  of 
the  incident  radiation  must  be  explicitly  provided. 

(4)  Choose  a  value  of  Teff  from  step  (2-3)  by  numerical  interpolation  or 
from  the  graph,  that  corresponds  to  the  incident  intensity  IQ .  This 
value  of  Tej£,  and  the  wavelength  of  the  incident  radiation, 
are  input  to  the  computer. 

(5)  The  computer  program  computes  and  plots  the  radiated  intensity  or 
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2 

power  spectra,  tha  intensity  of  the  incident  radiation  (watts/cm  ), 
and  the  analogous  blackbody  (Q^g-1 )  intensity  or  power  spectra. 

In  the  case  of  illumination  over  a  wavelength  spectrum  within  which  Qa^s(^) 
undergoes  substantial  variation,  Eq.  (2*16)  would  be  used  to  compute  vs 

jQabs(*)w(*)dX  over  the  incident  wavelength  spectrum.  The  latter  quantity 
is  the  average  intensity  of  the  incident  radiation,  and  would  be  substituted 
for  IQ  in  step  (3)  above. 

The  accuracy  of  the  numerical  Integration  routines  of  (3)  above  is  easily 
checked  by  letting  Qabg-1  (blackbody)  and  comparing  the  numerically 
evaluated  integral  with  the  analytic  result  of  the  Stephan- Boltzmann  Law: 
JVbb<A,T)dWT4. 


2 . 6  Results 


In  order  to  carty  out  the  computations  outlined  above,  we  must  choose  a 
theoretical  model  for  calculating  the  wavelength  dependent  absorption 
efficiency  factor  Qabs(*)‘  for  simplicity,  we  will  choose  for  this 
discussion  the  Extended  Quasistatic  (E.Q.S.)  model  ,  which  is  described  in 
Section  I.  This  model,  although  not  as  exact  as  the  Variational  technique 
(discussed  t.lso  in  Section  I) ,  has  the  advantages  of  being  analytic  and  of 
utilizing  sh^rt  computer  run  times.  In  the  case  of  engineering  applications, 
the  E.Q.S.  model  would  be  first  used  to  optimize  the  various  parameters. 
Having  accomplished  this,  one  would  then  use  the  Variational  technique  for  the 
final  computations. 

We  first  note  that  (E  parallel  to  particle  axis)  in  the  parameter  regime  where 
the  absorption  cross  section  is  maximized  (L«“<1  or  L'e"<l,  as  the  case  may 
be) ,  the  expression  for  the  absorption  efficiency  factor  reduces  to 

Q.bs  -»  o'  (U"  ->  0).  (2-18) 

Note  that:  (1)  in  this  regime,  the  absorption  efficiency  factor  is 
proportional  to  the  particle  radius,  a,  and  is  independent  of  wavelength. 
This  fact  can  be  utilized  in  the  tailoring  of  the  spectral  emissivity  of  a 
particle  cloud  e.g.  graphite  or  iron.  (2)  The  ratio  aabs/Qa|jS  *s  t^e 
projected  area  of  the  absorbing  particle,  with  respect  to  the  direction  of 
the  incident  radiation.  The  results  in  Fig.(l)  were  obtained  by  assuming  that 
the  incident  E-vector  is  parallel  to  the  particle  axis.  The  implications  of 
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this  will  U  discussed  in  a  subsequent  subsection.  (3)  For  the  assumed 
cylindrical  particlea,  tha  absorbad  power  is  proportional  to  A^-2al,  while 
the  radiated  power  is  proportional  to  A-2wal. 

2.6.1  Graphite. 

♦ 

Our  first  examples  of  snail  particulate  radiation  involve  particles  composed 
of  graphite.  The  first  example  will  use  particles  having  very  small  radii  of 

100  Angstroms.  Such  particles  can  be  produced  by,  for  example,  chemical  vapor 

2  5 

deposition  .  The  conductivity  of  these  is  estimated  to  be  10  mho/m.  The 

estimates  of  mean  free  time  and  mean  free  path  of  1.35X10' 15  sec.  and 

4X10’*^  m,  respectively,  were  provided  to  N.  Pedersen  by  I.  Spain,  Colorado 

State  University. 

Figure  1  is  a  log- log  plot  of  the  absorption  efficiency  factor  vs  wavelength 
for  the  graphite  particles  discussed  above.  The  particle  length  was  .'.hcien  to 
be  10  microns.  Note  that  the  coordinates  each  span  five  orders  of  magnitude. 
The  wavelength  range  is  from  O.lp  to  1cm.  Note  that  the  absorption  is  a 
constant  over  the  range  from  about  5p  to  about  5  millimeters.  Note  also  that, 
for  this  particular  particle,  <1  at  all  wavelengths  shown. 

Table  2-1  and  Fig. 2-2  were  computed  and  plotted  ,  as  previously  outlined,  by 
means  of  Inserting  a  temperature  Teff  in  Eq.  (2-17)  and  solving  for  the 
incident  intensity  IQ.  Note  that  the  wavelength  of  the  incident  beam  was 
chosen  to  be  10.6  microns.  This  computation,  involving  numerical 
integration  of  the  LHS  of  Eq.  (17),  was  carried  out  for  each  value  of  T  rr. 

Bit 

shown  in  Table  1.  Also  shown  in  the  Table  and  in  Fig. 2-2  are  the 
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corresponding  values  of  IQ  for  an  equivalent  blackbody  (Qa^g  "1  at  all 
wavelengths).  Since,  as  stated,  QajJS<l>  we  wou^  ®*pect  that  Te^.^  for 
this  particle  would  be  greater  for  a  given  value  of  IQ  than  ?eff  for  the 
equivalent  blackbody:  The  efficiency  for  absorption  is  greater  at  10.6  microns 
than  it  is  at  the  shorter  wavelengths.  Therefore,  the  particle  must  be  hotter 
(relative  to  a  blackbody)  in  order  to  radiate  the  same  power  that  it  absorbs. 

From  the  above,  it  is  clear  that  it  is  only  the  shape  of  the  Qa^s(^)  vs  * 
curve,  and  not  its  magnitude,  that  determines  the  curve  of  Te££  vs 
Io(Aq).  Table  1  and  Fig. 2-2  are  not  strictly  required  for  these 
computations,  but  are  provided  to  lend  insight  into  the  physics  of  the 
interactions .  . 


In  Fig. 2-3  are  shown  log-log  plots  of  spectral  emittance  (in  the  usual  units 
o 

of  watts/cm  -micron)  vs  wavelength  in  microns.  The  solid  line  represents 
radiation  from  the  graphite  particle,  and  the  dashed  line  is  the  usual  Planck 
blackbody  spectral  radiation  curve.  The  input  parameters  are  2500 
deg.K,  and  the  wavelength  of  the  incident  radiation,  Ao-10.6/i,  The  incident 
intensity,  Iq(Ao)-2600  watts/cm  ,  is  automatically  computed. 


The  blackbody  radiates  at  least  an  order  of  magnitude  more  efficiently  than 
the  particle.  This  is  because  Qa^s  is  on  the  order  of  0.1  or  less  in  the 
radiation  wavelength  regime  shown.  The  blackbody  spectrum  is  skewed  somewhat 
towards  the  shorter  wavelengths.  This  is  because  the  absorption  efficiency 
decreases  with  decreasing  wavelength  for  wavelengths  smaller  than  -5/i. 
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LOG (ABSORPTION  EFFICIENCY  FACTOR)  vs  LOG (WAVELENGTH,  ml  crons) 
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Figure  2-1.  Absorption  efficiency  factor  vs  wavelength  for  a  100 
Angstrom  radius  graphite  fiber. 
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Figure  2-2. 


Plots  of  T  w  vs  I  from  Table  2-1. 
ett  o 


A  -10.6  microns, 
o 


( 


58 


L  »«d»<  t lie  >■  1 0 . 6 
K»d  I  u*«  1E-8«i 
L«ngt  h"lE-5» 

Bulk  cond*lE3»!'0"i» 
m.  f .  p.  '<E-1Jii 
T*u«t.33E-13*«c 
6  Jul  1987 


hlcr omi 


Io 

<««U»/c»a2) 


lo<bl  «ekbodv*> 


T«ff 
(d*g.  k. ) 


o 


c 


• . 43E-0 1  - 

1.83E-01  - 

2.34E-01  * - 

2.97E-01  - 

3.76E-01  - 

4.77E-01  - 

6.03E-01  - 

7.62E-01  - 

9.62E-01  - 

1 . 21E+00  - 

1 . 53E+00  - 

1.92E+0O  - 

2.42E+O0  - 

3.03000  - 

3.88E+00  - 

4. 73£*0O  - 

3.93E+0B  - 

7.39E  +  00 - 

9. 1 9E+00  - 

1.14E+01  - 

1  *  4  1  Ei-01  - 

1 . 75E+0 1  - 

2.16E+01  - 

2. 66E+01  - 

3.27E+01  - 

4.01E+01  - 

4.91E+01  - 

6.00E+01  - 

7.32E*01  . 

8.9OE+01  - 

I . 08E+02  - 

1.31E+82  - . 

1 . 38E*02  . 

1.90E+02  - 

2. 29E*02  - 

2.73E+02  - 

3.29E+02  - 

3.93E*02  - 

4.69E+02  - 

5.39E+02  - 

6.64E+02  - 

7.87E+02  - 

9 . 3 1 E*02  - 

1. 1 0E*03  - 

1 . 30E+03  - 

1.33E+03  - 

1.79E+03  - 

2.  1 1E+03 - 

2.47E*C3 - 

2.88E+03  - 


1.44E-01  - 

1.84E-01  - 

2.3SE-01  - 

3.01E-81  - 

3.84E-01  - 

4.90E-G1  - 

6.26E-01  - 

8.00E-01  - 

1.02E+00  - 

1.30E+00  - 

1.67E+00  - 

2.I3E+00  - 

2.72E+00  - 

3.47E+00  - 

4.43E+00  - 

3.66E+P0  - 

7.23E+00  - 

9.23E+00  - 

1  .  18E  +  01  - 

1.30E+01  - 

1.92E+01  - — 

2 • 43E  +01  - 

3.13E+01  - 

4 . 00E+01  - 

3.11E+C1  - 

6.33E+01  - 

8.34E+81  - 

1.06E+02  - 

1.36E+02  - - - 

1 . 74E  *02 . . 

2.22E+02  - 

2. 83E+02  . 

3.62E+02  - 

4.62E+02  - 

3.90E+02  - 

7.53E+02  - 

9.62E+82  - 

1.23E+03  - 

1.37  0  03  - 

2.00E+03  - 

2.56E+03  - 

3.27E+03  - 

4.  17E  +  03 - 

5.33E+83  - 

€  80EI33  - 

8.69E+03  - 

1 . 1 1 E  +  04 - 

1 . 42E  *04 - - 

1 . 81 E»04  - 

2.31E+04  - 


300 
319 
339 
360 
383 
407 
433 
460 
489 
520 
353 
388 
£25 
664 
706 
731 
798 
848 
90» 
959 
1019 
1083 
1 151 
1221 
1301 
1383 
1470 
1363 
1662 
1767 
1878 
1996 
2122 
2236 
2398 
2349 
2718 
2881 
3063 
3236 
3461 
3679 
391  1 
1138 
4420 
4698 
4995 
3389 
5644 
6003 


I 

c  \ 


Table  2-1. 


Radiating  temperature  vs  incident  intensity.  Thin  graphite  fiber. 
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Figure  2- 


Figure  2- 


.  Radiant  emittance  of  the  100  Angstrom  radius  graphite  fiber. 
Effective  radiating  temperature  -  2500  deg.  K. 


4.  Radiant  emittance  of  the  10C  Angstrom  radius  graphite  fiber. 
Effective  radiating  temperature  -  1500  deg.  K. 
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We  next  lower  the  temperature  to  1500  deg.K  (Fig. 4).  The  required 

incident  intensity  at  AQ-10.6p  is  correspondingly  reduced  from  2600  to  52 
o 

watts/cm  .  The  radiation  spectra  of  the  particle  and  the  equivalent 
blackbody  are  significantly  reduced  in  magnitude  and  are  shifted  towards  the 
longer  wavelengths.  The  curves  are  more  nearly  equal  in  magnitude  because 
Qabs(*)  is  ®ore  nearly  equal  to  unity  for  these  radiating  wavelengths  than 
in  che  case  of  Teff"1500  de8- 


Figures  2-5,  2-6,  and  2-7  are  similar  to  Figures  2-3  and  2-4,  except  that  the 
values  of  (and  the  corresponding  are  lowered  to  1000,  500 
and  300  deg.K,  respectively.  The  ratio  of  particle  emittance  to  that  for  the 


blackbody  is  essentially  the  same  for  T@f£-  500  deg.K  and  300  deg.K.  This  is 
because  Q  is  essentially  constant  over  the  radiation  spectra  for  both 


values  of 
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Figure  2-5.  Radiant  emittance  of  the  100  Angstrom  radius  graphite  fiber. 
Effective  radiating  temperature  -  1000  deg.  K. 


Figure  2-6.  Radiant  emittance  of  the  100  Angstrom  radius  graphite  fiber. 
Effective  radiating  temperature  -  500  deg.  R 


It  is  next  instructive  to  see  what  happens  to  the  thermal  radiation  spectrum 
when  Qabg(*)  exhibits  more  structure  chan  that  shown  in  Fig. 2-1.  For  this 
purpose,  we  choose  a  graphite  particle  of  the  more  conventional  variety, 
having  a  length  of  1  millimeter  and  a  radius  of  1  micron.  At  the  short 
wavelengths,  the  correspondingly  large  values  of  (ka)  give  rise  to  the 
resonance  shown  in  Fig. 2-8.  Note  also  that  the  enhancement  of  Qa^g.  as 
predicted  by  Eq.(18)  is  seen  at  A— 3  millimeters. 

Figure  2-9  (IQ  vs  Tef£>  is  analogous  to  Fig. 2-2,  but  shows, as  expected,  a 
much  larger  deviation  from  the  blackbody  curve  of  Fig. 2-2. 

The  plots  of  radiation  spectra  for  this  particle,  with  (and  therefore 

IQ)  as  the  incremented  parameter  are  shown  in  Fig.s  2-10  through  2-13.  In 
the  case  of  this  particle,  it  is  easy  to  see  how  the  form  of  Qa^s(^) 
reflected  in  the  radiation  spectrum.  It  is  particularly  interesting  to 
observe  that,  at  the  highest  temperature,  i.e.  Te^f-2500  (Fig. 2-10),  the 
particle  radiates  as  though  it  were  much  cooler,  as  compared  with  the 
equivalent  blackbody.  At  the  lower  temperatures,  e.g.  Teff-500  (Fig. 9),  the 
radiation  spectrum  (neglecting  the  cusp)  has  essentially  the  same  shape  as  a 
blackbody  at  the  same  radiating  temperature. 
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Figure  2-10.  Radiant  emittance  of  the  1  micron  radius  graphite  fiber. 
Effective  radiating  temperature  -  2500  deg.  K. 
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Figure  2-11.  Radiant  emittance  of  the  1  micron  radius  graphite  fiber. 
Effective  radiating  temperature  -  1000  deg.  K. 
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Figure  2-12.  Radiant  emittance  of  the  1  micron  radius  graphite  fiber 
Effective  radiating  temperature  -  500  deg.  X. 


Figure  2-13.  Radiant  emittance  of  the  1  micron  radius  graphite  fiber 
Effective  radiating  temperature  -  300  deg.  K. 


2.6.2  Icon 


Also  of  interest  are  very  fine  filaments  of  Iron.  We  therefore  include  herein 
an  example  of  the  radiative  properties  of  such  a  particle  having  a  radius  cf 
100  Angstroms.  For  purpoaea  of  Illustration  of  the  potential  for  broadband 
screening  (via  absorption) ,  ve  have  selected  a  filament  length  of  100  microns . 

Figure  2-14  displays  the  absorption  efficiency  as  ?  function  of  wavelength. 
Note  that  the  curve  is  quite  similar  to  that  .Jot  100  Angstrom  radius  graphite 
(see  Fig.  2  -1),  except  that  the  maximum  value  of  Q  »^s  is  100X  that  for 
graphite.  This  is  because  the  electrical  conductivity  of  the  iron  is  taken  to 
be  <rc-10^  mho/m,  while  that  for  graphite  was  assumed  to  be  105  mho/m. 

In  Fig.  2-15,  we  see  that  the  plot  of  T#££  vs  IQ  is  essentially  the  same  as 
that  for  the  100  Angatrom  radius  graphite.  This  is  because,  as  mentioned 
previously,  the  shapes  of  the  Qabs(*)  curv®s  (hut  not  their  magnitudes) 
are  the  same. 

In  Fig.s  2-16,  2-17,  2-18,aftd  2-19  the  values  of  Tfi££  a*®  1500  (the  maximum 

permissible),  1000,  and  500  deg.K,  respectively.  The  following  comments  are 

of  interest:  <l)The  shapes  of  the  radiation  spectra  are  essentially  the 

same  as  for  the  thin  graphite  example  (Fig.s  2-4  through  2-7).  (2),  The 

required  incident  intensity  IQ  is  nearly  the  same  in  the  iron  vs  graphite 

o 

cases  (38  vs  52  watts/cm  at  T#££-1500  deg.K)).  The  reason  for  this  is 
apparent  from  inspection  of  Equation  (2-17).  (3)  The  magnitude  cf  the 
emittance  is  greater  by  a  factor  of  -100X  for  the  iron  vs  the  graphite. 
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Figure2-14.  Absorption  efficiency  factor  vs  wavelength  for  e  100 
radius  iron  fiber. 
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Figure  2-15.  Plot  of  Tg££  vs  Io  for  a  100  Angstrom  Iron  fiber.  Aq-10.6/). 
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Figure  2-16.  Rediant  emittance  of  ths  100  Angstrom  radius  iron  fiber 
Effective  radiating  temperature  -  1500  deg.  K. 


LOG C EM JTTAMCE ,  cron )  v*  LOG C WAVELENGTH,  micron*) 


Effective  Radiating 

Temp."  1800  deg.K 

Lamdr.C  me  1  dant  )>18. 

S  micron* 

IngrVtftiit  tntana'ty- 

1 .  1E+B1  ei*tt*/'cm-2 

/  \ 

-Blaekbody  Spectrum 

/  \  \ 

\  \ 

\  \ 

\  \ 

\  \ 

\  \ 

\  \ 

\  \ 

\  \ 

\  \ 
\  \ 
\ 

\  PARAMETERS  Cmka) 

\ 

\  R*dtu*"lE-8  m 

\ 

\  Length* 1E—4  m 

\ 

\  Cond.“lE7  mho/Tn 

\ 

\  T*u“4E-l5  aec 

\ 

\ 

\  m.f.p.«3E-ll  m 

‘ulm 


Figure  2-17.  Radiant  emittance  of  the  100  Angstrom  radius  iron  fiber 
Effective  radiating  temperature  -  1000  deg.  K. 


Figure  2-18,  Radiant  enictance  of  the  100  Angstrom  radius  iron  fiber. 
Effective  radiating  temperature  -  500  deg.  K. 


Figure  2-19.  Radiant  emittance  of  the  100  Angstrom  radius  iron  fiber. 
Effective  radiating  temperature  -  300  deg.  K. 
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2 . 7  Applications 


In  this  subsection,  we  consider  some  interesting  applications  of  the  foregoing 
treatment.  These  include  various  radiative  and  absorptive  effects  associated 
with  electromagnetic  radiation  incident  upon  aggregates  of  microscopic 
absorbing  particles.  We  will  also  consider  the  conditions  required  for 
vaporization  and/or  melting  of  such  particles.  The  environment  is  assumed  to 
be  exoatmospheric .  Before  addressing  these  subjects,  however,  it  is  necessary 
that  we  discuss  the  geometrical  considerations  which  arise  when  the  absorbing 
particles  are  nonspherical  and  randomly  oriented. 

2.7.1  Geometrical  Considerations 

As  already  mentioned,  the  foregoing  analyses  were  carried  out  under  the  tacit 

assumption  that  the  particle  axis  is  aligned  parallel  to  the  electric  field 

vector  E  of  the  incident  wave.  This  of  course  presents  no  problem  for 
o 

spherically  symmetric  particles.  However,  in  the  case  of  nonspherical 

particles,  we  are  interested  in  the  absorption  efficiency  factor,  as  a 

function  of  the  angle  of  incidence,  as  well  as  the  polarization  angle.  We  are 

not  permitted  to  take  the  average  value  of  Instead,  the 

absorbed  power  must  be  calculated  using  the  quantify 
2  2 

cos  (a  )  Qab  s  (  *  ,  f  £  )  "Co  s  (  Q  )  °  at,  c  (  *  >  *  i  )  /  2  a  1  s  i  nS  ^  ,  where  0^-angle 
of  incidence  and  a-polarization  angle.  The  thermal  radiation  on  the  other 
hand,  must  be  calculated  as  before,  using  the  total  particle  area  2iral. 

Each  value  of  6^  will  therefore  give  rise  to  a  different  radiation  spectrum. 


72 


As  can  be  seen,  the  precise  calculation  of  thermal  radiation  from  a  small 
ensemble  of  randomly  oriented  nonspherical  particles  involves  a  somewhat 
different  treatment  than  that  of  Section  2.3.  The  following  accounts  for 
arbitrary  particle  oriental  on  in  the  rase  of  our  E.Q.S.  theory: 

We  firs.  a  that  the  absorption  cross  section  is  proportional  to  the  square 
cf  the  component  of  the  incident  E- field  that  is  projected  on  the  particle 
axis.  Thus,  including  polarization  angle  a  and  the  angle  of  incidence  6^, 
^he  absorption  cross  section  becomes 

<Tabs^A,0t,ei^  "  £rabsi:A*0>’r/2)cos2(Q)sin2^i)-  (2‘19) 

The  absorption  efficiency  factor  is  inversely  proportional  to  the  area  of  the 
particle,  projected  in  the  plane  perpendicular  to  the  incident  k  vector. 
Thus , 

(WA’a’®i)  -  aabs(A>Q’ V/2alsin(*i>-  (2'20) 

Combining  Eq . ' s  (2-19)  and  (2-20),  the  absorption  efficiency  factor  for 
arbitrary  particle  orientation  (E.Q.S.  theory)  is 

QabsvA'a>Si)  “  %bs(A'0'*/2)cos2(a)sin(*i)/2al-  (2-21) 

In  the  case  of  monochromatic  incident  radiation,  the  absorbed  power  is 

F  .  -I  (A  )A  Q  ,  (A  ,a ,9.) 
abs  o'  o'  p^abs'  o  ’  i 


“  Io^A0^p<7abs^  Ao  ’  0  ',r/2)cos2  (a)  sin(0  i)/2al , 


(2-22) 


and,  in  the  present  nomenclature,  the  radiated  power  is 


«0 

rad  “  AJQabs ( ‘ X  ’ 0 • */2 ) ^5 (e (hc/XkT 


2*c2h 


-dX. 


(2-23) 


eff'  -1) 


o 

For  cylindrical  particles,  A/Ap-2iral/2al-ir.  For  monochromatic  Illumination, 
an  expression,  analogous  to  Eq.  (2-17),  results  from  equating  Eq.’s  (2-22)  and 


(2-23): 


Qabs(Ao,0,w/2)cos2(a)sin^i) 


"i 


2wc2h 


I  mbs  A5(e'hc/AkTeff>  -1) 


dX. 


(2-24) 


We  see  that  Eq.  (2-24)  Is  transcendental  with  respect  to  the  variables 
a  and  0^. 

The  problem  can  be  solved  numerically  as  follows: 

(1)  Assume  values  of  the  intensity  and  wavelength  of  the  Incident  beam.  i.e. 

assume  that  we  know  I  (X  )  . 

ov  o' 

(2)  Successively  increment  the  angles  a  and  0^  to  include  all  possible 
orientations  of  the  particle  with  respect  to  the  incident  E-vector. 

(3)  For  each  of  the  above  angular  increments,  calculate  L.H.S.  of  Eq .  (2-24). 


(4)  For  each  above  angular  increment,  solve  Eq.  (2-24)  (by  numerical  iterative 


means)  to  compute  the  correct  value  of  Teff 

(5)  Having  obtained  Te££,  compute  the  radiation  spectrum  for  that  specific 
combination  of  a  and 

(6)  Repeat  steps  (1)  through  (S)  for  each  angular  Incrementation.  Add  the 
spectral  intensities  (vatts/cn>2)  according  to  wavelength  for  all  angles. 

The  above  procedure,  using  appropriate  angular  weighting  (sin(0^))  and 
normalization  for  the  numerical  integration  of  (6),  will  yield  the 
orientation-averaged  thermal  radiation  spectrum  from  our  cylindrical 
particles.  With  modification,  the  more  exact  Variational  technique  for 
calculating  Qaljs  can  be  used  in  the  above  method. 


2.7.2  Aggregates  of  Absorbing  Particles 
Infernally  Heated  Spherical  and  Cylindrical  Clouds 

We  now  consider  the  case  of  a  cylindrical  cloud  of  particles.  Assume  that  the 
source  of  illumination  is  uniformly  distributed  along  the  cylinder  axis,  and 
that,  for  all  wavelengths  of  interest,  the  product  7(A)R»1,  where  R-radius  of 
the  cylinder,  and  y-absorption  coefficient  of  the  cloud  given  by 

y(A)  -  n*abs(A)  m’1,  (2-25) 

.3 

where  n-number  density  of  absorbing  particles  'm  ) . 
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In  this  csss,  the  input  radiation  per  unit  axial  length  Is  simply 

Pin(tot)/L,  where  P^n(tot)-total  Input  power  and  L-length  ot  the  cylinder. 

The  radiated  power  oust  equal  the  input  power,  and  is  emitted  over  a  surface, 

or  "skin"  whose  depth  is  approximately  7“*.  The  effective  </alue  of  the 

emlssivity  of  the  cylindrical  cloud  is  Just  that  of  a  blackbody.  This  is 

because  of  the  following  argument:  (1)  Vithin  the  cloud  (r«(R-y  A)  and 

r»7**),  the  particles  are  essentially  in  thermal  equilibrium  with  their 

neighbors  and  the  radiation  spectrum  in  such  regions  is  essentially  blackbody 

radiation.  (2)  The  effective  number  of  particles  that  are  radiating  at  the 

surface  at  any  wavelength  is  N-nA  7(A)'^,  where  Ap  -total  cloud  area. 

c  c 

Thus,  the  larger  QajjS (^ )  is  at  some  wavelength,  the  smaller  the  number  of 
radiating  particles  that  are  radiating  from  the  "surface"  at  that  wavelength. 
(3)  Therefore,  the  radiation  from  the  surface  of  the  cloud  will 
be  essantinlly  blackbody  radiation  and  the  effective  emlssivity  of  the  cloud 
will  be  unity. 

The  radiating  temperature  TQ^£  of  the  cloud  can  therefore  be  found  from  the 
following  equation. 


If  the  cloud  were  spherical,  and  the  illuminating  source  were  placed  at  the 


canter  of  the  sphere,  then  the  sane  procedure  can  be  used  by  replacing  the 

o 

factor  2*RL  in  Eq.  (2-26)  by  4*R  ,  where  in  this  case  R-sphere  radius. 

This  equation  demonstrates  an  interesting  (and  potentially  useful) 
circumstance  that  an  optically  thick  (yR»l)  absorbing  cloud  of  particles  will 
radiate  essentially  as  a  blackbody  whose  surface  temperature,  and  therefore 
its  radiation  spectrum,  can  be  controlled  by  its  surface  area. 

Cylindrical  Cloud;  External  Illumination 

We  next  consider  the  case  of  an  axially  illuminated  cylindrical  cloud  of 
absorbing  particles.  This  problem  only  has  meaning  if  7R«1.  If,  in  addition, 
7l>l,  then  the  incident  radiation  within  the  cloud  is  given  by 

W'>  -  V1*.  <2-2'> 

where  z  is  along  the  cylinder  axis.  The  particles  will  radiate  with 
emisslvlty  Q^s(X)W^(A,z) .  Thus,  the  radiation  spectra  of  all  surface 
elements  of  length  dz  would  be  added  according  to  wavelength.  As  long  as  the 
above  criteria  on  7R  and  7L  obtain,  the  overall  spectrum  from  the  cylindrical 
cloud  will  depend  only  upon  Qa^g(A) . 

If  the  7(A)L  were  «1,  then  all  particles  would  radiate  in  accordance  with  Eq. 
(2-17),  usirg  I^nc-IQ(Ao) ,  and  e»Q  ^  (A)  ..  *n  tne  case  of  axial  solar 
illumination,  Eq .  (16)  is  applicable  here  for  computing  the  radiation 
spectrum. 
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Slab  Geometry;  External  Illumination 


The  case  of  external  illumination  of  a  slab  composed  of  absorbing  particles 
can  be  treated  simply  (see  Ref.  24,  p.513)  using  the  diffusion  approximation. 
This  is  appropriate  when  -)rr»l  and  when  a  .  »  a  ,  where  r-slab 
thickness.  Ue  further  assume  that  Qa^s(^)  is  fairly  uniform  over  the 
wavelength  ranges  of  interest.  The  problem  is  to  find  the  temperatures  at  the 
illuminated  front  surface  T  ^(1)  and  at  the  back  surface  Teff(2)>  which 
is  assumed  to  be  not  illuminated. 


Siegel  and  Howell  provide  a  simple  calculational  procedure  for  the  computation 
of  the  above  temperatures.  Their  results  can  easily  be  used  in  conjunction 
with  our  computational  methods  to  determine  the  radiation  spectrum  emanating 
from  both  surfaces. 

We  first  note  that,  if  it-*,  no  thermal  radiation  will  occur  at  surface  2  and 

the  incident  intensity  I^nc  must  therefore  be  radiated  from  surface  1,  which 

will  radiate  as  a  blackbody.  Therefore,  the  S tef an- Boltzmann  law 

(I, _ j-I  -oT^)  to  calculate  the  maximum  temperature  at  surface  1.  Let 

raa  o 

this  be  Te££(max) . 

We  now  let  yr  be  fairly  large  but  finite.  The  temperatures  at  the  two 

24 

surfaces  are  then  given  by 

Teff(1)  "  Teff(max) [((l/2)+(3/4)7r)/(l+(3/4)7r)]1/4,  (2-28) 

and 


Te f f ( 2 )  -  Teff(max) t(l/2)/(l+(3/4)70 ] 1/4 . 


(2-29) 
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Equation  (2-2)  can  now  be  used  to  calculate  the  blackbody  spectra  emanating 
from  the  two  surfaces.  As  long  as  the  slab  absorbs  virtually  all  of  the 
incident  radiation,  the  spectral  distribution  of  this  will  not  sensibly  effect 
the  spectra  emanating  from  the  two  surfaces. 

Incident  Electromagnetic  Pulse 

There  is  interest  in  the  interaction  between  an  absorbing  particle  cloud  and  a 
strong  electromagnetic  pulse  (BMP) .  Assume  that  the  pulse  shape  is  square  and 
has  a  duration  of  10  microseconds.  This  means  that  the  predominant  power  in 
the  pulse  is  in  the  frequency  range  of  from  100  kHz  to  ~500kHz  (300m>A«>60m)  . 
In  order  for  the  particles  to  be  strongly  absorbing  at  these  and  higher 
frequencies,  the  electrical  conductivity  should  not  be  excessive,  and  the 
aspect  ratio  of  the  particle  must  be  quite  large.  For  example,  if  we  were  to 
select  the  100  Angstrom  radius  particle  of  Fig.  (2-1)  for  this  purpose,  its 
length  would  (from  the  E.Q.S.  theory)  have  to  be  2:2  mm  and  the  filament  would 
have  an  aspect  ratio  of  1Q5,  if  this  could  be  accomplished,  and  if  the 
filaments  were  not  aglomerated,  this  pulse  would  be  strongly  absorbed  by  a 
cloud  of  such  particles.  The  scattering  would  be  negligible.  The 
temperature  would,  because  of  the  small  radius,  be  uniform  throughout  the 
volume  of  each  particle. 
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Vaporization  and  Melting 

Under  conditions  of  extremely  intense  incident  illumination,  particles 
composed  of  graphite  will  vaporize.  Taking  the  vaporization  temperature  of 
graphite  to  be  4200  deg.C  -4473  deg.K,  and  using  our  first  example  (e.g.  Fig. 
(2-1))  of  a  thin  absorbing  graphite  fiber,  we  find  that  the  intensity  of  an 
incident  beam,  of  wavelength  Ao-10.6/i,  is  1300  watts/cm?. 

Sufficiently  thin  iron  particles  will  melt  under -much  lower  illumination 
levels  than  give  rise  to  the  vaporization  of  the  above  graphite  fiber.  This 
is  because  (1)  the  melting  point  of  iron  is  lower  (1808  deg.K  vs  4473  deg.K), 
and  (2)  a  sufficiently  thin  filament  of  iron  has  a  much  larger  absorption 
efficiency  factor.  If  wa  take  the  iron  particle  discussed  in  Section  2.6.2  as 
our  example,  and  Ao-10.6/i,  we  find  that  an  incident  intensity  of  only  65 
watts/cm?  is  required  to  melt  the  subject  particle  in  a  space  environment. 

We  believe  that,  upon  melting,  the  iron  filament  would  (due  to  surface 
tension)  form  a  spherical  droplet.  If  this  is  correct,  the  spherical 
particle  would  no  longer  be  highly  absorbing  and  would  cool  and  solidify.  It 
is  therefore  not  at  all  likely  that  any  vaporization  would  occur  in  the  case 
of  the  iron  particle,  except  in  the  case  of  extremely  high  intensity  levels. 

Ionization 

For  the  vaporized  graphite  fiber,  we  would  have  (we  believe)  a  gas  of  neutral 
carbon  atoms.  Assuming  this  to  be  the  case,  a  field  strength  on  the  order  of 
108  to  10$  volts/cm  would  be  required  to  give  rise  to  field  ionization  of  the 


carbon  atoms. 


If  either  of  the  above  particles  were  present  in  a  gaseous  environment  (e.g. 
the  atmosphere)  ,  field  ionization  of  neutral  molecules  would  occur  at  very 
much  lower  incident  field  strengths.  This  is  due  to  the  fact  that  very 
substantial  enhancements  of  the  incident  E- field  will  occur  very  near  the  tips 
of  the  filaments.  This  interesting  topic  has  been  the  subject  of  considerable 
in-house  work  at  Panametrics.  We  might  also  note  that,  if  conductive  highly 
elongated  particles  were  embedded  in  a  solid  dielectric,  we  predict  that  the 
incident  intensity  level  required  to  initiate  dielectric  breakdown  would  be 
vary  substantially  reduced. 


Required  Total  Mass 

Also  of  interest  is  the  mass  required  to  accomplish  a  given  (d£)  reduction  in 
intensity  of  an  incident  beam.  This  calculation  is  straightforward  in  the 
case  of  a  slab  geometry:  The  total  mass,  Mtot>  in  the  aggregate  is 

Htot  -  Acn,V'  <2'W) 

.  3 

where  A  -illuminated  cloud  area,  n-number  density  of  particles  (m  ), 
c 

r-depth  of  the  cloud  (m) ,  Vp-p»rticle  volume,  and  p-mass  density  of  the 
material  comprising  the  particle.  We  note  that  (1)  n-7/eabg,  aiM*  (2)  the 
one  way  (dB)  attenuation  (absorption)  is  given  by  (dB)-(2 . 3/10)yr .  Combining 
these  and  Eq.  (2-28),  we  have  the  expression  for  the  cloud  mass  per  unit  area 
required  to  cause  a  given  (dB)  reduction  of  an  incident  beam: 

Mtot/Ac  "  <2-3/10)(dB)pVp/%bs(Ao).  (2-31) 
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2.8  Coawnti 


JJ ,,|NTU!ll  UJJI..I.   -J  ,  ....  . -  ■.  ,J 


The  foregoing  discussions  relating  incident  elect  romagnetic  energy, 
microscopic  particle  parameters,  and  the  resulting  thermal  radiation  were 
undertaken  to  provide  a  general  understanding  of  the  Important  processes  as 
they  relate  to  the  physical  parameters  of  the  particles.  The  specific 
problems  discussed  were  chosen  for  mathematical  simplicity  and  ease  of 
physical  understanding.  The  standard  texts  (e.g.  Ref's.  24  and  25)  on 
Radiation  Theory  can  bo  used  in  conjunction  with  the  work  of  this  Section  in 
the  solution  of  more  complicated  problems  -  and  there  are  many. 

The  thermal  emlasivlty  of  small  particles  has  been  discussed  by  Kattawar  and 
26 

Eisner  who  also,  we  note  with  pleasure,  discuss  the  quantum  mechanical 

27 

aspects  of  the  emlsslvlty  for  small  particles.  Pluchino  has  computed 

emissivlties  (equal  to  Qgj,,  and  less  than  unity)  of  small  layered  spherical 

particles.  Our  initial  objective  herein  has  been  to  address  the  question  of 

26 

«(*)"QaKa(A)»l.  We  also  note,  however,  that  Bohren  and  Huffman  have 
addressed  this  question  and  are  in  agreement  with  our  conclusions  that  this 
circumstance  is  perfectly  reasonable.  They  also  provide  some  Interesting 
historical  aspects  of  the  subject. 

Other  objectives  of  this  Section  have  been  to  compute  the  appropriate 
radiation  spectra  for  various  specific  particles  of  interest,  and  to  determine 
means  of  computing  the  radiative  properties  of  aggregates  of  broadband  highly 
absorbing  particles  under  conditions  of  intense  illumination. 
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Affttiii  A 


SCATTBAIM  BI  CVKTIB  COMMTCTITS  VIBK1S 
Pxil iainary  Iictij 

All  the  ABitclotl  results  obtained  by  Pima*  tri  c  t  to  the 
pniiat  tin  for  loattiriai  and  absorption  by  conductive  fibers 
have  dealt  vitk  straight  fibers.  The  question  naturally  arises 
then  as  to  .the  effeets  of  curvature  of  the  fibers  on  the 
scattering  and  absorption  efficiency. 

A  search  of  the  literature  reveals  very  little  work  on 
curved  wires;  all  of  that,  with  two  exceptions,  involving 
perfectly  oonduoting  wires.  The  curved,  perfectly  conducting 
wire  was  apparently  first  considered  by  Aharoni  in  1946  . His 
equations  were  applied  to  circular  loop  and  spiral  antennas  by 
Msl.ll  In  1936  Kosyosaj ian  considered  ba  c  k- a  oa  t  te  r  ing  from 
perfectly  conducting  oircular  loops. I1  The  two  exceptions  to'the 
perfectly  oonduoting  case  are  the  work  of  Philipsoa,  who 
considered  lossless  dieleotric  rings, and  Aoquists,  who 
considered  wavy  cylinders. 30  in  both  of  these  latter  oases, 
however,  the  scatterer  was  taken  to  be  only  a  perturbation  on  its 
surroundings,  so  that  a  full  integral  equation  approach  was  not 
required. 

Ve  have  derived  the  integral  equation  for  curved  fibers, 
having  finite  conductivity,  fros  first  principles.  The  usual 
thin-wire  analysis  invariably  assuaes  that  the  electric  field  can 
be  expressed  in  teras  of  a  current  filaaent  concentrated  on  the 
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•  sit  of  th*  fiber,  tad  this  it  physically  sonewhat  unsatisfactory 
especially  when,  it  in  the  present  cate,  we  not  t  work  with  both 
eleotrio  and  nagnetio  fields  at  tka  cnrfaoe  of  the  fiber.  Ve  use 
instead  an  approaok  bated  on  Hnygan't  principle, 31  vbiok  states 
tkat  fields  generated  by  tke  tangential  oonponents  of  S  and  H 
(distribntad  along  tke  snrfaoe  of  tke  fiber)  must  precisely 
oanoel  tke  axial  oonponents  of  tke  incident  electric  field  along 
tke  fiber  axis. 

This  results  in  an  equation  involving  integrals  of  the  two 
unknown  functions  E  and  H  along  tke  fiber.  Taking  tke  thin—  wire 
Unit  where  fiber  radius  is  very  snail  oonpared  to  incident 
wavelength,  tke  seoond  of  these  integrals  is  fairly  straight¬ 
forward,  and  for  good  eonductors  is  interpretable  as  tke  field 
due  to  a  distribution  of  surfaoe  currents.  Tke  first  integral 
behaves  differently,  however.  Tke  kernel  reduces  to  a  delta- 
function,  resulting  in  a  tarn  in  E  at  the  field  point  of 
evaluation  of  tke  integral  equation.  Snrfaoe  values  of  E  and  H 
can  then  be  related  by  a  surface  inpedanoe  conoept  to  finally 
give  a  pure  integral  equation  for  the  ourrent. 


Tkt  Iit«|til  Ifutioi 


When  an  incidant  •  1  ec  t  r  on  agn  •  ti  o  field  E*nc  illuninates  * 
body  in  free  spaoe,  the  resulting  fields  are  related  by  Huygens' 
principle,  which  states  rigoronsly  that^l 

Bin0(r)  -  ( 1/  4it )  yxjdo'ki'  z  E+ (  r  '  )  g(kH) 

1E(r)  outside  (/.-D 
0  inside  . 

The  left-hand  side  (LHS)  of  this  equation  consists  of  the  suai  of 
the  eleotric  fields  due  to  the  incident  wave  and  surface 
distribution  a  over  the  body  of  Magnetic  and  eleotric  dipoles, 
respectively. 

Here 


g(kR)  -  (l/kR)e*«,  R«  I  r  -  r '  I  , 


( A-2 ) 


where  t_'  and  r  are  the  souroe  point  and  field  poi.it, 
reape  otively . 

Equation  (A-l)  states  that  the  B  field  is  given  by  the  LHS 
for  all  field  points  outside  the  surface.  On  the  other  hand,  for 
all  field  points  within  the  surface  the  LHS  vanishes  identically, 
i.e.  the  surface  field  distributions  sust  precisely  cancel  the 
incident  wave.  This  latter  stateaent  is  sosetisei  known  as  the 
extinction  theorem,  or  the  extended  boundary  condition. 
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The  extinction  theorem  it  applied  to  the  curved  fiber,  shown 
in  Fig.  A-J,  as  follows.  Let  a  **  fiber  radius,  and  p(s)  “  radius 
of  curvature  as  a  function  of  position  a  along  the  axis  of  the 
fiber.  We  assume  that 

•  / Pain  <  <  1  ( A-3 ) 

a/b  <  <  1  ( A-4 ) 

ka  <  <  1  ,  ( A-5 ) 


i.e.,  the  fiber  radius  is  mucL  leas  than  the  minimum  radius  of 
curvature,  the  fiber  half-length  b,  and  free  space  wavelength  k  = 
2n/X,  respectively.  Ve  also  assume  the  fiber  to  have  moderate  to 
large  conductivity,  so  that  axial  currents  will  be  induced  and 
guided  along  the  fiber. 

Now  requiring  that  the  axi  al.  oomocuent  of  Eq.  (.4-1)  vanish 
along  the  fiber  axis  gives 

( 1/  4n )  a  •  yx  ^x  J*  do '  in'  x  H+  (  r.'  )  g  (kR) 

+  (l/4«)  s  •  ^x^dcr'  kn  '  x  E+  ( r '  )  g(kR)  =  s  •  E*no(s)  .  (A-6) 


Note  that  this  equation  is  still  exact,  although  we  have  only 
used  a  portion  of  the  information  available.  The  curl  operators 
may  be  taken  under  the  integral  sign,  because  the  field  point 
need  never  approach  the  fiber  scrfaise. 
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it  ai  i  a<i  ari  aet  anaitaMUni 


fe  suppose  the  magnetic  field  on  the  surface  to  be  purely 
azimuthal,  so  that 


2>  <*'  )  H(  s  '  )  , 

«'  *  (t')  -  H ( s * )  . 


(A-7) 


Using  the  identity  enrl  curl  ■  grad  div  -  div  grad,  the  first 
integrand  of  Bq.  (A-6)  takes  the  torn 

s  •  v*  V*  ft’  *  H+  ( )  g(kR)  -  s  •  (7  V  —  V*  V )  s '  B(e')  g(kR) 


-  B(s')  s  *7V-  s'  g( kk)  +  k2(s*a')  H(s')  g(kR) 

-  H(s')  (s.V)  (S'*V)  g(kR)  +  k*  (*  •  s')  H( s ' )  g(kR) 

-  -H(s')  [a^/g.g.t  -  k2  (s  •  s')  g ( kR)  ,  ( A-d) 


where  in  the  last  step  we  hare  used  the  fornal  notation  s  • V  - 
d/ds  and  s  1  ‘  V  ■  -i/d  s'  for  direotional  derivatives.  The  minus 
sign  arii.  s  in  the  latter  case  because  the  primary  variable  has 
the  form 


2  *  r  {  s )  -  r  (  s  '  ) 

For  the  second  integral  of  Eq.  (A-6),  the  electric  field  on 
the  surface  is  assumed  to  be  purely  axial,  i.e. 


s'  B( s  1  )  . 


(A- 9) 


Nov 


•  •  Vx  n '  x  B+(r' )  g(kR)  -  - «  •  (n'  xE+)  x  Vg 

■  ( n '  xE+)  •  ( »  x  vg ) 

■  E(t’)  (n'  x*')  •  (i  x  vg) 

■  -  B  (  a  ' )  ( *  •  s')  (  n '  •  vg) 

-  -E( a ' )  (s  •  l')  (£'  •  R)  k  g ' 

“  -B(s')  (s  •  ;•)  ( a / R) k g '  .  ( A- 10 ) 

Hare  in  the  fourth  atop  ve  used  the  identity 

( a_  x  b)  •  ( c.  x  d_)  “  ( a  •  c )  ( •  d)  -  (a*  d)  ( b.  •  c.) 

then  noted  that  the  scalar  product  n’  •  s  vanishes  identically 
under  the  azimuthal  portion  of  the  surface  integration.  In  the 
last  step  above  g'(kR)  ■  dg/d(kR)  and  ve  have  assumed  that 

n'  •  &  =  a/R  .  ( A-l 1 ) 

Note  froai  Pig*  A-l  that  this  equality  only  holds  when  I  s- s '  I  << 
p,  i.e.,  for  points  sufficiently  close  together  along  the  axis 
that  the  curvature  cf  the  fiber  has  not  come  into  play.  No 
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Approximation  it  involved  here,  however*  because  the  integral  is 
entirely  negligible  otherwise,  at  we  will  now  see. 

Using  Eq.  (A-10) ,  the  second  term  of  Eq.  (A-6)  becomes 


-U2/4n) 


E(s* ) (•'•') 


(a/E)  g'UR) 


-  - ( k/2 ) 


Us) 2 


E( a  1 ) (s- s') 


(1/kR)  g’UR) 


+• 

*  (l/2)a2  B( a ) j  4«'[(»-i')2  +  a2]"3'2  -  Eta)  .  (A-12) 


From  the  first  step  one  notes  that  because  ka  <<  1  only  that 
portion  of  the  integral  of  order  (1/ka)*  will  contribute  to  the 
final  result.  Noting  that 


(1/kR)  g'UR)  »  -(1/kR)3  =  -  ( 1  /  k ) 3  Ks-s’)2  +  a2]“3/2  , 


for  kR  <<  1,  the  remaining  steps  of  Eq.  (A-12)  are  straight¬ 
forward. 
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Patting  the  results  of  Bqs.  (A-8)  and  (A-12)  back  in  Bq.  (A' 
6)  and  oarrying  out  th  a  aiiauthal  integration,  ona  finds  that 


( i a/2 ) 


ds'H(s') [d^/ds 


8  a ' -k1 ( s • s ' ) ] 


g(a, s'  ) 


B(  s ) 


-  -i  .Eiae  (a)  .  (A-13) 

Va  oan  axpraaa  thia  raanlt  aa  an  integral  aquation  for  total  1  lna 
ourtant  1(a)  by  writing 

1(a)  -  2nai  H(a)  .  (A-14) 

and  introducing  a  surface  inpadanca  par  unit  length  given  byl* 


Z  -  -  ( e>|i/2nakf )  J0  ( kf  a)  /  J  0 '  (  kf  a)  .  (A-15) 

Hera  kf  -  (»2|i*  +  i u  p  o) 2  is  the  oonplaz  propagation  constant 

within  the  fiber,  J0  is  the  Bessel  function  of  the  first  kind, 
and  ws  have  assumed  that  both  E  and  H  within  the  fiber  vary  inch 
aore  rapidly  in  the  radial  than  the  axial  direction.  This 
assuaption  is  consistent  with  the  reqoireaent  of  aoderate  to 
large  conductivity.  Note,  however*  that  if  one  were  to  represent 
I',)  Ks  a  Fourier  expansion,  then  eventually  the  axial  variations 
of  such  terns  would  doninate.  with  the  result  that  k£  in  Eq.  (A- 
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IS)  weald  have  to  be  Modified. 1$  That  ia,  th •  lurfaoi  inpedance 
baooaaa  dependent  oa  tba  cata  of  axial  variation  of  currant  whan 
that  rata  ia  larga. 

Tha  tangential  B  field  can  now  be  expressed  aa 

E ( a )  -  ZI ( a )  .  ( A- 1 6 ) 

and  using  thia  raanlt,  along  with  Bq.  (A-14),  one  finally  obtaina 


Note  that  for  atraight  fibera  a*a*  »  1  and  thia  aquation  reduce* 
to  tha  nanal  fornnla.5  4iao  for  curved,  perfectly  conducting 
wirea  Z~  *0  and  one  again  find*  tha  accepted  foranla.l^ 


Alftoiiaitioa  £•(  Sfioial  Pikers 


For  the  general  csss,  ss  described  by  tbs  integral  Bq .  (A- 
17)  ,  it  i  <*  clear  that  detailed  nuaerical  computations  era 
required  in  order  to  obtain  any  erplioit  results.  If  the  fiber 
axis  baa  radius  of  curvature  large  coapared  to  wavelength,  or  is 
aade  rp  of  a  xigsag  aeries  of  atraigbt  line  segments,  however, 
then  the  absorption  cross-section  is  readily  appr oxiaa ted  using 
earlier  results. 

First,  oonsiier  the  case  where  the  fiber  axis  curves  only 
slowly  and  it  relatively  long,  ooapared  with  wavelength.  Then  we 
can  use  the  infinite  oylinder  result  to  obtain  the  absorption 
cross-section  per  unit  length  o,/*,  where  this  ratio  is  given  by 
Bq.  (1-20).  Note  that  o,  will  be  a  function  of  position  a  along 
the  curved  fiber,  in  that  ea  depends  on  the  angle  60  foraed  by 
the  incident  B  vector  and  the  looal  tangent  to  the  fiber  axis. 
The  total  absorption  cross-section  ^*bs  then  given  by 

♦  b 

£,bs  "  I  da  to,  (•)/*]  .  (A-18) 

-b 

It  is  now  straightforward  to  obtain  results  for  toroidal  or 
C-shaped  fibers,  or  other  configurations  aeeting  the  slowly- 
crrviug  limitation,  by  numerical  integration  of  Bq.  (A-18). 

For  a  fiber  aade  up  of  xigsag  line  segaents,  it  is  sort 
appropriate  to  eaply  the  quasi-static  appr oxia a ti on  to  the 
absorption  cross-section  o,b,  »*  given  by  Eq.  (1-29).  We  have 
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already  mi  by  coipitiioi  with  other  eoapntations  that  this 
(«mli  is  qaite  looititi  for  oil  fiber  lengths  provided  tho 
aodified  depolarisation  faetor  L*  of  Bq.  (1-31)  io  md  for  k0£  > 


Nov  for  i  fibor  of  N  segaeats,  kit  it|  individual  lengths  tu, 
n  ■  1,  1.  ...  ,  N,  tha  total  absorption  eroaa-aootioa  is 

given  siaply  by 


^aba  *  £  «aba  <*a>  •l*2  «n  • 


( A-19) 


skin  tho  factor  sis^g,  (angle  between  incident  B  veotor  and  nth 
segaent  azia)  is  iaoladed  beoanae  the  original  Bq.  (1-29)  was 
speeifloally  for  breadaido  inoldoaoo. 


It  is  netful  to  note  that  for  both  Bqa.  (A-lg)  and  (1-19) 
i<  or oaa-seetioa  is  a  liaoav  ana  of  tho  or oie-ssetioni  of 


iadividnal  segaents.  B< 


>h  sack 


pr eei sely 


as  a  straight  fiber  nnde r  orientation  averaging,  one  oonolndoe 


that  the  orientation-averaged  oross-aeotio 


a  for  the  on *ve< 


fiber  will  be  jest  equal  to  that  of  a  straight  fiber  of  the  seas 
total  length. 


-jinanmanAr«£rtJu-uiruvutf\MV>AefNMJinAruwuwu 


A|fii4is  1 


TAMBT  OBSCVBATIOM 

la  this  notion,  taking  iato  ooa  »i  da  r  a  t  ioa  botk  tka  diffuse 
cioad  aeattariag  propartias  and  tka  kaaa  extinction  properties, 
oaa  requires  tka  optiaua  eoaditioaa  foe  targat  obaoaratioa  froa  a 
aaaa  efficiency  standpoint,  ovat  a  wavelsagtk  raaga  wkere  tka 
tkeory  ia  believed  to  ba  aoat  aocarata.  It  ia  appropriate  to 
aaploy  tka  quasiatatie  approilaatioa,  ia  wkiek  partiolaa  ara 
aaaaaad  to  ba  aaall  ia  soapariioa  vitk  aaaaor  wcvelengtk. 

Coaaidar  tka  aitaatioa  wkea  a  targat  of  (radar  or  optioal) 
or oaa-aact ioa  oj  ia  partially  eoaoaalad  by  a  eload  of  partiolaa, 
tka  oload  having  aa  iaookaraat  cross-section  Ec,  Because  of 
tkia  iaookaraaoa,  tka  oroas-seotioas  ara  additive  aad  oaa  kaa 

Observed  oroaa-aaotioa  ■  +  oy  ,  ( B- 1 ) 

wkara  y  aad  x  ara  tka  daoay  eoaitaat  aad  tkiokaaaa  of  tka  eload, 
respectively. 

Wo  aow  atiua  tkat  tke  aaaaor  oaaaot  detect  tka  target  if 

2C  i  I «T  e-2yx  (B-2) 

wkara  C  is  a  figure  of  aerit,  e.g.  for  K  ■  0.1  tke  aigaal  froa 

tke  target  would  ba  10  dB  down  in  tke  "  not  at"  of  tke  elond 


raters 


Tki  diffesa  letttiriii  fro*  tbs  cloud  kn  aarliar  k«ci 


ooapetsd  to  b«32»S3 


E 


c 


(B-3) 


*kin  Vp  is  si  agl  s-pa  r  ti  c  1  s  volui,  tk«  rsdaesd  coadaotivity  is 
givsa  by  os '  “  kQs  -  Vp0/«0  oc  ( ia  tins  o f  coadeetiv  ity ,  or  tbs 
imsgiaaxy  part  of  tbo  relative  disisotrio  eoastsat).  sad  Ac  is 
tbs  geometrical  orosr-sictioa  of  tbs  cload.  or  tbs  ssasor  bsaa 
eroas- sse t  i  oa  at  tbs  oloid,  if  tbs  lattsr  shoald  bs  saalisr. 
Usiag  this  resalt,  aad  takiag  tbs  sqaality  ia  Bq .  (B-2)  tbsa 

giVSS 


*  1  *1  S“2*T 


(B-4) 


Tbs  dsoay  ocastaat  for  tbs  oload  is  gives  by 


T  ■  B  • 


(B-3) 


ia  tsras  of  tbs  aaabsr  dsasity  of  partiolss  aad  tbs  orieata  tioa- 
avs  rigsd  sxtiactioa  or o s s- sse t i oa  psr  partiols.  Assaaiag  that 
absorption  sffsots  doaiaats,  oas  bas 

*o  V 

®*xt  6  5  sbs  **  - ® r -  o  Oq  V»/  3  (  B-6 ) 

311  +  (Ls  ri 

wksrs  ia  tbs  asxt-to-last  stsp  tbs  qaasiststic  spproxiaati oa  was 
saploysd,  aad  ia  tbs  last  stsp  we  aotsd  that  optimum  sbsorptioa 
will  ooear  whoa 
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L«"  -  4  ( •/£)  2  l  aU/a)  -1]  <<  1 


(B-7) 


(L  is  tka  depolarisiag  factor).  Not*  site  that  tka  total  sail  M 


of  pi 


*  ■  P  V,  h  A.  . 


(B-8) 


Baplojiii  Bqs.  (B-J).  (B-6)  aid  (B-t)  ia  Bq.  (B-4),  and  taking 

tka  logaritka  of  tka  raaalt  detsraiaes  tka  required  aaai  to  ba 


( B-9) 


Tka  quaatitiaa  K.  aad  Ae  ara  praacribad  by  tka  logistics  of 

tka  task  at  kaad.  Tkua.  to  aiaiaixe  N  oaa  first  aust  saal.  a 
aatorial  vitk  anal laat  possibla  value  of  tka  ratio  p/oc'.  Ia 
additioa,  it  is  desirable  to  obtaia  as  large  a  value  of  tka 
product  kQ2  Vp  a Q '  aa  possibla.  ia  order  to  aiaiaixe  the  loga- 
ritkaic  tara  ia  Bq.  (B-9). 


For  coacr ataaass.  suppose  tkat 


Li"  -  0.1  . 


k„*  -  0.1 


(  B-l 0) 


Froa  Bq.  ( E-7 )  we  thsr  find  that 


(B-ll) 

particle  radiaa,  once  the 
Typical  particle  design 

a)  microwave  ration 
y  -  3 .x  10“2  a 
l  -  4.8  x  10-4  B 
oc  »  5  xlO®  aho/a 


a  -  2.2 

x  10“8  a 

inf  cared 

reaion 

y  ■  3  x 

10-5  . 

l  -  4.8 

X  10-7  b 

«c  “  1  x 

104  aho/a 

a  “  1 .0 

x  10-8  a 

0.01  (I/a) 

4  UnU/a)  -1) 


This  equation  can  be  nsed  to  determine 
length  and  conductivity  axe  known, 
parameters  can  now  be  listed: 
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n m  mkji  r\a  r.*  iul  muum  rxWlXJUI  UfurjUMAKlUilUtll  rfUUi  TLk  iLrflttfTLM 


Apftiilz  G 


TAILOKING  OF  PABTICLB  PABAMBTBB8  FOB  BPBCIFIC  APPLICATIONS 

(The  discussion  in  this  Appendix  vss  done  during  1  985  ss 
psrt  of  the  effort  of  the  first  jreit  of  the  Progress  The 
theoretical  and  computational  techniques,  as  described  in  Section 
1  of  the  present  report,  have  been  si  gni  f  ioantly  improved  since 
that  time.  However,  this  Appendix  is  included  in  order  to 
demonstrate  the  approach  we  have  devised  for  this  application.) 


It  is  the  purpose  of  this  section  to  demonstrate  the  tech¬ 
niques  by  which  the  parameters  (length,  radius,  and  conductivity) 
of  thin  oondactive  fibers  can  be  adjusted  ao  that  a  cloud  of 
these  particles  will  have  selected  specified  spectral  charac¬ 
teristics.  These  characteristics  are:  (1)  large  absorption  and 

small  scattering,  (2)  large  scattering  and  low  absorption,  and 
(3)  transparency.  Two  frequencies  were  arbitrarily  chosen  for 
the  oases  to  be  analyzed.  These  are  fj  =  101®  Hz  and  fj  »  10*2 
Hz,  corresponding  to  wavelengths  of  3  cm  (microwave)  and  300 
microna  (infrared). 

The  analyses  to  be  discussed  will  be  based  on  our  "extended 
quasistatic  model",  the  basis  of  which  is  described  in  Section 
1.3.  A «  demonstrated  in  Section  1,  this  anslytical  model 
produces  results  which  are  in  surprisingly  good  agreement  with 
the  results  obtained  using  the  much  more  rigorous  variational 
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technique  of  Section  1.2.  To  be  acre  explicit,  the  queeistatic 
model  can  be  need  to  calonlete  the  absorption  crcee  aection  over 
a  very  vide  range  of  k0i,  including  k  Ql  >>  1.  The  scattering 
croae  aection  computations  resulting  fron  the  quaeiatatic  model 
are  valid  only  in  ^he  range  k0£  n. 


In  all  of  the  cases  analyzed,  the  results  cf  the  quasistatic 
calculations  are  directly  compared  with  the  corresponding  results 
of  the  variational  technique.  Ve  find  these  comparisons  to  be 
quite  remarkable. 


CotitUitiv*  Iftttioti 


la  addition  to  tho  equations  of  Sections  1.2,  1.4  and  1.5. 
the  following  aquations  arc  utilized  in  tha  foregoing  analyses: 

Orients t ion- averaged  absorption  cross-section 
(see  Eq.  1-29): 

®  abs  ■  °abs  (C-l> 

Orientation-averaged  scattering  c r oa s- a e c t i on  (see 
Eq.  1-30): 

®aca  “  (j)  ®aoa  (C-2) 

High  freqnenoy  depolarising  factor  (see  Eq.  1-31): 

L‘  -  4  (-J-)  Un  -1]  (k0£  <  n)  ( C— 3 ) 

The  inclusion  of  the  above  equation  poraits  extension  of  the 
salculaticn  of  absorption  cross  seotion  well  into  the  kQ£  >>  l 
regine.  Due  to  this,  ve  refer  to  the  present  theory  as  the 
"extended  quasistatic  theory." 

In  this  treetaent,  ve  are  using  the  Drude  (free  electron) 
model  for  the  coaplex  dieleotric  constant  (Section  1.4)  and  the 
Fuchs  model  for  the  dependence  of  conductivity  on  partiole 
dimensions,  and  ve  have  purposely  chosen  frequencies  and  wj 

such  that  possible  anomalies  duo  to  these  do  not  appear.  Also, 
ior  purposes  of  analysis,  ve  can  sake  the  simplification  e  * 
is".  The  computer  program,  however,  does  not  utilize  this 
siapl if i ca  t icu. 
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Absorb  st  3  so,  loflsst  st  300  aisroms 

Tin*  criteria  lead  to  tba  following  nathaaatical 
atataaanta : 


Le  t 

Juju.  .  .01 

(C-5) 

•  p  (2> 

(La"  (1)  )2  <<  1 

( C-6 ) 

Let 

La"  (1)  -  0.1 

(C-7) 

Equation  (C-7)  snsnzaa  that  the  absorption  oroaa  section 
will  be  at  its  aaziaua  valve  st  frequency  ,  but  will  be  reduced 
by  a  large  footor  at  trequeney  *2  ”  2*  z  1022  Hx .  Also.  we  know 
that  the  scattering  cross  section  will  be  at  least  as  great  at  wj 
as  at  .  This  will  causo  otCa(2) /«abs^2^  to  oe  large,  which  is 
what  we  seek  at  wj • 

Froa  Eqs.  (C-5),  (C-7),  (C-l),  and  (C-2),  wo  obtain  (with 

the  appropriate  definitive  equations  of  Section  1.3)  equations 
tor  the  conductivity  «c,  radius  a  and  length  h 
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! 


1/3  .  1/2 

o  (  .03)  (o) _ 

< til)17*  (10  Un  (f)  -11  }1/6(«c)1‘r5 


t 


2a  10{Uu  (j)  -1] 


( C-9 ) 


(C-10) 


*h*a  Bqs.  (C-8)  and  (C-9)  are  equivalent,  i.e.,  v a  on  choose  a 
eonduotivity  and  aolva  for  a.  or  choose  a  radius,  a*  and  solve 
for  tho  conductivity  oc.  Va  will  do  tka  forair.  Knowing  botk  a 
and  oc ,  wa  tkan  calculate  tka  langtk  from  (C-10).  Nota  that  tka 
tarn  in  aqnaro  braokata  in  (C-8),  (C-9),  and  (C-10)  ia  vary 
insansitiva  to  larga  variationa  in  (t/ a).  Tkia  ia  aapaoially 
troa  wka n  it  ia  raiaad  to  fraotional  powara  in  (C-8)  and  (C-9). 
Tharafora,  wa  can  coma  quita  oloaa  to  tka  daairad  raanlta  if  wa 
siaply  lat  [  ]  *  2.5.  Tkia  peraits  diraot  estimation  of  a  and  t, 
given  a  praaalaotad  value  of 


Table  I 


Cond 

m 

1  .00E+04 

Radiua  “ 

8.68E-06 

Langtk 

at 

3  .7  IB-03 

k0(l) t/2  - 

3 . 8  8E-01 

Cond 

m 

1 . 00E+05 

Radiua  “ 

2  . 7  5  E-  0  6 

Langtk 

m 

3  .71E-0  3 

k0(lU/2  - 

3 . 88E-01 

Cond 

m 

1  . 00E+06 

Radius  ■ 

8.68B-07 

Langtk 

m 

3 .7  IE-03 

k0(l)*/2  - 

3 . 8  8E-01 

Cond 

m 

1 .00B+07 

Radius  “ 

2  . 7  5  E-  0  7 

Length 

m 

3  .71B-03 

k0(lU/2  - 

3  .8  8E-01 

Cond 

a c 

1 . 00  E+0  8 

Radiua  * 

8.68E-08 

Length 

“ 

3  .7 IE- 03 

k0(l)£/2  - 

3  . 8  8  E-  0 1 

Cond 

m 

1 .00E+09 

Radius  “ 

2.75  E-08 

Langtk 

«* 

3  .7  IE-03 

k0(l)£/2  - 

3  . 8  8E-01 
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The  above  pt  ooidnn  leads  to  sets  of  peraissible  parameter 
values,  as  oxeaplified  in  Table  I.  Substitution  of  the  various 
parameter  aats  into  the  ooapater  prograa  yields  graphs  suoh  as 
shown  in  Fig.  C-l .  Sinoe  these  graphs  are  identical  over  the 
wavelength  range  300  aiorons  to  3  oa,  only  one  is  shown.  This 
figure  is  representative  of  a  fairly  highly  conducting  aetal 
(bulk  oonduotivity  -  10?  aho/a) .  The  relaxation  tiae  t  is  that 
of  ooppar. 

Proa  Fig.  C-l,  we  see  that  (1)  indeed,  the  absorption  orosa 
section  peaks  at  very  nearly  the  wavelength  X*  »  3  oa,  (2)  that, 
in  the  oaae  of  absorption,  the  extended  quasistatic  theory  is  in 
excellent  agreeaent  with  the  variational  theory  for  wavelengths 
greater  than  about  30  aiorons,  and  (3)  the  scattering  as 
calculated  froa  both  theories  ia  in  excellent  agreeaent  when 

M  *  1  • 


We  further  conclude  that  the  aatheaatioal  "design" 
prooedure  deaoustrated  in  this  sub-seotion,  although  not 
analytically  exaot,  provides  good  paraaetrio  sets  for  the 
solution  of  the  stated  problea,  since  the  objectives  are  net  at 
the  iwo  specified  wavelengths. 


Figure  C-l .  An  example  of  predoainant  absorption  at  X  ■  3  ca 

and  pradoainant  scattering  at  300  aicrona.  v  “ 
variational  taoniqua;  a.q.s.  *  Extended 
Qnaaiatatic  Theory. 
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Akittk  it  I  limb  at  100  ■& 


Vi  kivi  quit*  different  otitirii  for  this  problsa  ss 
ooapsrsd  with  the  preceding  problsa.  Bars,  ws  wish  to  have  s 
vary  highly  absorbing  olond  which  has  low  soattiriig  oross 
•action. 


In  this  case* 
given  by 


tha  clond  scattering  cross  saotion 


is 


£ 


o 


,4C 


9sca 
*  «.bs 


) 


$ 


(C-ll) 


whore  A0  ■  projected  area  of  the  olond  or  that  portion  thereof 
which  ooonpies  the  solid  angle  of  the  incident  bean. 

Obviously.  we  wish  to  aake  the  absorption  oross  section  of 
the  particle  nooh  larger  than  its  scattering  oross  section.  This 
anst  be  true  over  at  least  the  wavelength  range  f roa  300  aicrons 
to  3  ca. 

The  above  considerations  lead  to  the  following  aatheastioal 
criteria : 

(1)  In  order  to  have  aasinua  absorption  over  the  specified 
wavelength  range,  we  want  Ls"  (1)  <<  1  and  Ls"  (2)  << 

1.  This  aeans  that  the  psrticlos  shonld  be  very  thin 
and  (perhaps)  not  too  highly  conducting. 
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(2)  li>M,  aaiir  tko  ikevi  oaiitloai,  *•««  **  k0*  Vp*  «Q2 
•*<  -  Vj  «0.  tk«  ratio  (lig|/ltki)  it  proportional 

to  k02  Vj  «e.  fkoroforo,  ia  ordor  to  aaiataia  kigk 
eload  akaorptioa  aad  lev  eload  aoattotiag.  wo  witk  to 
koop  tko  prodaot  Vf  *0  low  (wo  kawo  ao  ooatrol  owor 
ke). 

A  littlo  kit  of  oapo riaoati ag  witk  tko  oatoadod  gaaaiatatie 
(MS)  ooapator  prograa  yiolda  appropriato  iota  of  paraaotora. 
Tko  paraaotcr  aot  oorroapoad lag  to  aiaiaal  total  aaaa  it  a  -  100 
Aagatroaa,  t  -  100  aicroai,  aad  «0  -  10*  ako/a  at  tko  kaTk 
ooadastiwi tp .  Tko  r»talta  of  aaiag  tkoao  paraaotora  ia  tko 
BOS  aad  wariatioaal  oodoa  aro  akowa  ia  Fig.  C-2 .  Oao  caa  trado 
off  *>  laxgor  r^diaa  for  a  lowor  eoalaotiwitp  aad  aekiowo  tiailar 
roauta.  kat  witk  a  aoaowkat  kigkor  ratio  of  (0««a/dBg9). 

% 

Noto  agaia  tkc  rraarkaklo  agrooaoat  kotwooa  tko  eoapatatioat 
katod  oa  two  two  ladogaadoat  tkoorioal 


Fi|tt •  C-2 .  Aa  ixiaplt  of  broadbaad  nbaorptioa  at  3  oa  aad  300 
aloroai.  y  -  variatioaal  toehaiqaa;  a.q.a.  - 
Bsttidid  Q  aaistatio  Theory. 
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va  ry  low  tbiotf tion 


lof laatiag  Mm 

Xt  it  mijt  to  <oil|i  fiititlti  kivlai 
m4  relatively  kith  aaottoriag.  Bs|m1m«»  koo  ahowa  we  that  if* 
for  maflo.  wo  ehooio  i  highly  uoadaotiag  metal  omoh  to  eoyyar 
or  alam  laws*  oad  ooloot  i  ttiiai  of  ooo  oitroa  or  grootor*  tho 
tortiolo  will  ho  ooooatiolly  aoayletaly  roflootiog  for  ot  leaat 
oil  wowolomgtht  ho low  90  aieroaa.  Vo  ooa  tailor  tho  aaottoriag 
oroaa  oootioa  to  hoooao  yroyortioaol  to  a*  at  wavcioagtha  aigaif~ 
ioaatly  ho low  that  for  whiah  k0£  *  l»  *hiah  of  eoaraa*  royroaoata 
laylalgh  aaottoriag  froa  '^erfaetly  ooadaetiag"  wiree.  Thao* 
tha  aggrogata  of  yartioloa  oaa  oaaily  ho  aada  to  ha  roflaatiag  at 
oaa  vavelaagtk  aad  aaaaatlally  traaayaraat  at  aaaa  (aigaifi- 
aaatly)  lowar  wavalaagth.  Chaff  aloada  hahawa  ia  thia  way. 
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It  tkia  out,  ve  «nt  tki  loittoriii  cron  tootiot  to  bo 
■•ok  lower  tkot  tko  obaorptioa  oton  sootioa  ot  nil  wovolot|tki 
aader  ootiilocttloi.  Tko  ro(tifoootto  Ito  to  tkio  otd  tko 
troaopareaej  oriterioa  ot  oj  tto  tioea  below: 

(1)  Noko  Vp  0Q  00  oaoll  oo  possible,  ouiol  otoot  witk  tko 
otker  criteria. 

(2)  Hoko  Lo  "  -  1  ot  m2. 

())  Moko  k 0t  <  1  ot  M2* 

Appl loo  tioo  of  tkooo  oritorlo  load*  to  tko  poroaetero  liotod 
lo  Pigoro  C-3 .  Noto  tkot,  ot  1  ■  )  oo.  tko  oboorptioo  it  tkroo 
orders  of  aagaitade  lower  tkoa  tkot  ot  X  -  300  oleroot.  Tko 
ooottoriat  tkroogkoot  tko  mi  to  m 2  opootroa  io  atok  lover  tkoa 
tko  oboorptioo. 


Ill 


log(Crota  SaetAinlt  vol) 


PLOTS  Of  ELECTROMAGNETIC  CROSS  SECTIONS  VS  WAVELENGTH 


Fi|«n  C-3  .  As  mapl*  of  trinpitiaoy  it  3  em  ltd  pr«dostii*t 
•biorptioa  at  300  aicroaa. 
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mmiMHMmuiiu  IUIUIUI 


CtMiiti  IcUtivi  to  Paxtiele  Paraaeter  Tailoring 


I 

h 


(1)  Our  extended  quasictatic  treituent  provides  s  very 
useful  tool  in  the  design  of  particles  for  specified 
applications.  In  sots  instances,  the  EQ S  theory  does 
not  yield  accurate  quantitative  data  and  should  be  used 
primarily  as  s  first  step  to  be  followed  by  the  full 
oosputation  using  our  Variational  aethod. 

(2)  The  calculation  of  s  oa  tterlns  cross  sections  using  the 
EQS  theory  is  limited  to  the  range  kc£/2  <_  1.  so  the 
variational  technique  aust  be  eaployed  beyond  that 
poi nt . 

(3)  We  have  included  backscatter  cross  section  computations 
in  Figures  C-l,  C-2  and  C-3  for  those  applications  that 
require  this. 

(4)  Within  the  bounds  of  present  technology,  it  is  possible 
to  create  particle  parameters  that  could  provide  useful 
spectral  characteristics  over  a  wide  wavelength  range. 
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